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Introduction 


This  report  is  composed  of  a  series  of  individual  technical  summaries  which  span  several  fun¬ 
damental  research  projects  associated  with  traveling  wave  charge  transport  technology  deve’  ip- 
ment.  A  detailed  description  of  the  GaAs  traveling  wave  charge  transport  device  and  its  potential 
capabilities  may  be  found  in  Reference  4.  This  device  is  of  interest  for  transversal  filtering  and 
other  signal  processing  applications  due  to  the  precise,  high  speed  charge  transfer  and  distributed 
charge  sensing  functions  which  can  be  accomplished  in  the  device  over  a  very  broad  bandwidth. 

The  projects  described  in  this  report  fall  into  three  basic  areas:  device  physics  and  .modeling, 
fabrication  technology  research,  and  GaAs  surface  acoustic  wave  (SAW)  research.  An  overview  of 
the  fundamental  technology  development  tasks  is  presented  via  a  brief  description  of  the  indivi¬ 
dual  projects. 

Device  Physics  and  Modeling 
1  >  raveling  W  ■  Charge  Injection  Research 

The  complicated  theoretical  problem  of  charge  injection  into  a  buried  channel  traveling  wave 
charge  transfer  device  is  addressed.  The  problem  is  analyzed  by  formulating  an  adequate  physical 
model  for  the  injection  process  and  solving  this  model  to  obtain  theoretical  predictions  and  perfor¬ 
mance  trends  for  injection  parameters.  This  research  involves  extremely  intensive  numerical  com¬ 
putation. 

2)  Mon-Destructive  Sensing  Research 

A  theoretical  model  is  formulated  to  describe  charge  detection  by  capacitively  coupled  sense 
electrodes  on  the  surface  of  the  traveling  wave  charge  transport  device.  The  detection  characteris¬ 
tics  peculiar  to  this  mode  of  signal  detection  are  studied. 
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3)  Transport  Research 

Various  field  problems  are  studied  and  analyzed  to  assess  the  impact  of  non-uniform  surface 
structures  in  the  charge  transport  device.  The  result  of  this  work  may  also  be  applied  to  the 
analysis  of  the  conventional  charge  coupled  device. 

Fabrication  Technology  Research 

4)  Low  Process  Temperature  Dielectric  Research 

This  project  involves  a  study  of  the  polyimide  dielectric  for  GaAs  1C  applications  where  low 
process  temperature  schedules  are  required  to  preserve  the  structural  and  electronic  quality  of  cir¬ 
cuit  features.  This  is  particularly  useful  in  applications  where  Schottky  barrier  diode  leakage 
degradation  due  to  heat  treatment  must  be  minimized. 

5)  Deep  Implantation  Research 

The  present  architectures  for  the  traveling  wave  charge  transfer  device  utilize  relatively  thick 
active  layers  for  electronic  channel  definition.  This  project  involves  the  development  of  photolitho- 
graphicallv  definable  masks  for  high  energy  ion  implantation  processing  of  thick  active  layers  and  a 
study  of  the  material  properties  resulting  from  high  energy  implantation. 

6)  Schottky  Barrier  Leakage  Research 

Charge  transfer  devices  require  low  values  of  channel  leakage  current  for  optimum  operation. 
This  project  involves  a  study  of  primary  factors  which  influence  reverse  bias  leakage  current  in 
GaAs  Schottky  barriers. 


GaAs  Surf  ace  Acoustic  Wave  Research 
7a)  GaAs  SAW  Waveguide  Research 

Acousti:  waveguides  provide  an  elegant  means  for  achieving  controlled  propagation  of  high 
power  acoustic  fields.  SAW  generation  and  propagation  in  GaAs  waveguide  structures  is  investi¬ 
gated  and  the  application  of  these  structures  in  the  traveling  wave  charge  transfer  device  is 
explored. 

7b)  GaAs  SAW  Propagation  Research 

The  achievement  of  transversal  filter  functions  in  the  traveling  wave  charge  transfer  device 
requires  SAW  propagation  through  large  arrays  of  surface  electrodes.  This  project  involves  a  study 
of  the  fundamental  properties  of  SAW  propagation  and  diffraction  in  both  bare  GaAs  and  GaAs 
with  surface  electrode  structures. 


1.  Traveling  Wave  Charge  Injection  Research 


Objective 


through  the  SAW  velocity  since  by  definition  current  is  charge  transported  per  unit  time.  There¬ 
fore  the  size  of  the  charge  packet  which  is  injected  into  the  transport  channel  for  a  specified  device 
bias,  signal  voltage,  doping  profile  and  SAW  potential  amplitude  is  the  essential  quantity  which  the 
model  will  provide.  To  calculate  this  quantity,  it  is  necessary  to  solve  for  the  mobile  carrier  distri¬ 
bution  of  the  injected  charge  using  the  coupied  set  of  nonlinear  differential  semiconductor  equa¬ 
tions.  See  appendix  A  for  details  of  the  solution  method. 

Due  to  the  nonlinearity  of  this  two-dimensional  problem,  it  is  necessary  to  utilize  a  numerical 
algorithm  to  solve  the  system  of  coupled  equations.  The  specific  algorithm  that  is  employed  for  the 
present  problem  was  developed  by  Herbert  Stone. ^  It  was  chosen  since  it  has  been  used  success¬ 
fully  by  other  authors  in  applications  involving  high  voltages  and  large  areasJ2,3^  which  are  typical 
conditions  found  in  the  traveling  wave  device. 

Theoretical  Results 

By  separating  the  injection  problem,  as  posed  in  Fig.  1,  into  several  simpler  cases,  basic 
processes  which  contribute  to  the  total  solution  are  made  evident.  This  also  allows  for  a 
verification  of  the  correctness  of  the  numerical  algorithm  without  the  complexities  of  the  complete 
problem.  The  algorithm  is  tested  with  two  cases: 

1)  An  input  section  with  no  SAW  potential. 

2)  A  transport  channel  far  removed  from  the  input  section  with  a  SAW  potential  impressed  on 

it. 

Figure  2(a)  shows  an  input  section  and  transport  channel  of  the  traveling  wave  device.  Bias 
voltages  are  applied  to  the  various  contacts  but  the  SAW  is  absent.  In  this  steady  state  condition, 
mobile  charge  surrounds  the  ohmic  input  contact  while  the  transport  channel  is  devoid  of  carriers. 
This  is  illustrated  in  Fig.  2(b)  which  is  a  calculated  contour  plot  of  the  charge  distribution.  The 
region  of  high  charge  density  surrounds  the  ohmic  contact  because  it  is  a  point  of  maximum  poten¬ 
tial  and  a  source  of  electrons.  Moving  away  from  the  contact,  the  carrier  density  decreases  and 
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Transition  Region 


eventually  approaches  zero  in  the  transport  channel.  The  transition  region  which  separates  high 
carrier  density  region  from  the  low  density  region  is  approximately  5  Debye  lengths  wide  and  is 
exhibited  on  the  plot  by  the  series  of  closely  spaced  equi-charge  density  curves.  They  are  labeled  as 
fractions  of  the  max:mum  density  which  exists  at  the  ohmic  contact.  The  transition  region  pro¬ 
duces  a  boundary  around  the  mobile  carriers  which  lacks  sharpness  due  to  carrier  diffusion. 

The  effect  of  a  SAW  potential  is  introduced  into  the  model  by  adding  a  sinusoidal  variation  to 
the  doping  density  profile  which  is  an  apj.  imate  rentes'''' UaUftci  of  t?ie  true  *  induced  polari¬ 
zation  charge  density.  The  confinement  of  mobile  charge  in  the  resulting  two-dimensional  potential 
well  is  clearly  exhibited  in  Fig  2(c).  This  carrier  density  plot  is  drawn  from  data  that  was  calcu¬ 
lated  at  a  location  far  removed  from  the  input  contact.  Again,  the  smeared  packet  boundary  is  evi 
dent.  The  problem  of  charge  bunching  in  a  SAW  induced  potential  well  has  been  previously 
solved.14*  The  solution  method  and  the  approximations  employed  differ  from  the  present  technique. 
However,  the  results  shown  here  agree  qualitatively  with  that  work.  Further  analysis  is  required 
to  correlate  the  two  solutions. 

The  analysis  of  the  injection  problem  involves  solving  for  the  carrier  distribution  which  exists 
at  the  input  section  when  the  SAW  potential  is  present.  This  is  essentially  a  joining  of  the  two  pre¬ 
vious  cases.  In  this  situation,  charge  is  trapped  in  the  SAW  potential  well  and  the  packet  breaks 
away  from  the  input  source  of  electrons  as  the  wave  progresses  into  the  transport  char,  tel.  The 
distance  which  the  SAW  travels  before  a  packet  splits  from  the  electron  source,  depends  on  the 
contact  voltage.  V,c.  It  is  convenient  to  measure  this  distance  from  the  center  of  the  input  contact 
where  the  SAW  is  assigned  a  phase  of  0  degrees.  Since  the  location  ot  split  off  is  not  known  before¬ 
hand.  an  iterative  procedure  is  performed  which  involves  solving  for  the  carrier  distribution  in  a 
series  of  steady  state  problems  of  increasing  SAW  phase.  The  phase  is  increased,  simulating  wave 
motion,  until  it  is  determined  that  a  charge  packet  has  split  off.  At  this  point,  the  amount  of  charge 
in  the  injected  packet  is  obtained  by  integrating  the  carrier  density  distribution  over  the  packet 


boundaries. 


This  analysis  is  performed  on  the  device  structure  shown  in  Fig.  1  for  an  8  /xm  wavelength 
SAW  in  a  4  /tn  thick  n-layer  that  is  uniformly  doped  with  a  donor  concentration  of 
Nd  =  1015/cm3.  A  bias  voltage  of  OV  is  applied  to  the  Schottky  electrodes  with  the  bottom  of  the 
n-layer  also  held  at  OV.  The  input  contact  is  biased  with  3V.  The  SAW.  now  with  its  true  y 
directed  variation,  is  specified  to  havt  a  1/2  V  zero  to  peak  potential  at  a  depth  where  the  channel 
electron  potential  resulting  from  the  depletion  charge,  the  bias  voltages  and  the  SAW  is  a 
minimum. 

The  computed  results  are  presented  in  Fig.  3.  Each  frame  shows  the  carrier  distribution  at 
instants  in  time  specified  by  fractions  of  T,  where  T  is  the  time  for  the  SAW  to  travel  one 
wavelength.  The  equi-charge  density  contour  lines  differ  by  1.4  X  1014/cm3  with  the  outer  most 
line  representing  a  concentration  of  7  X  1013/cm3.  At  t=0,  the  SAW  electron  potential  minimum  is 
located  at  the  center  of  the  input  contact.  The  resulting  carrier  distribution  is  shown  in  frame  1. 
In  frame  2.  the  SAW  has  propagated  1  Uj.i  to  the  right  and  pulls  the  mobile  carriers  with  it.  As 
one  potential  well  of  the  SAW  propagates  away  from  the  input,  a  preceding  well  is  advancing 
towards  the  input.  The  fields  of  this  well  attract  mobile  electrons  from  the  input  contact  and  con¬ 
sequently  it  will  fill  with  charge.  Each  of  these  two  wells  will  hold  the  same  amount  of  charge 
when  they  are  symmetrically  located  about  the  input  contact  as  shown  in  frame  3.  The  last  frame 
depicts  the  split  off  of  a  packet  from  the  electron  source.  For  this  bias  condition,  split  off  occurred 
after  a  time  of  .7 IT  =  1.98  ns. 

Since  ViC  varies  the  charge  load  of  the  injected  packet  and  consequently  the  channel  current,  a 
transconductance  curve  san  be  generated  by  repeating  the  above  analysis  for  different  values  of 
V|C.  This  was  performed  and  the  result  is  shown  in  Fig.  4  for  a  400  /xm  wide  channel.  The  param¬ 
eter  <f>m  is  the  minimum  electron  potential  du’  to  ionized  donors  in  a  layer  depleted  of  mobile  car¬ 
riers  as  calculated  using  a  1  dimensional  analysis.  If  VIC  is  less  than  <f>m,  a  conducting  channel 
would  exist  between  the  input  contact  and  output  contact  of  the  transport  device  in  the  absence  of 
a  SAW  field.  Therefore  <i>m  represents  the  minimum  input  contact  voltage  required  to  deplete  the 


n-layer.  By  subtracting  this  value  from  the  applied  contact  voltage  as  done  in  Fig.  4,  it  enables  the 
transconductance  of  devices  with  different  bias  and  doping  parameters  to  be  conveniently  com¬ 
pared 

Conclusion 

A  model  of  the  charge  injection  mechanism  in  a  GaAs  buried  channel  traveling  wave  charge 
transport  device  is  developed.  It  includes  the  effects  of  substrate  parameters  and  electrode 
geometry.  Results  are  generated  for  two  special  cases  which  are  used  to  test  the  numerical  algo¬ 
rithm.  The  potential  and  mobile  carrier  density  distributions  in  a  biased  transport  channel  are 
obtained  in  the  vicinity  of  the  injecting  contact  when  the  SAW  is  absent  and  at  a  location  removed 
from  the  injector  when  the  SAW  is  present.  The  process  by  which  a  charge  packet  is  injected  into 
the  transport  channel  is  demonstrated  and  a  characteristic  transconductance  curve  for  the  traveling 
wave  charge  transport  device  is  calculated. 


2.  Non-Destructive  Sensing  Research 
Objective 

The  objective  of  this  research  is  to  study  the  fundamental  properties  of  non-destructive  sens¬ 
ing  of  charge  in  a  traveling  wave  charge  transport  device  by  providing  a  theoretical  model  for  the 
small  signal  time  and  frequency  domain  responses  of  the  non-destructive  sense  structure  in  terms 
of  the  physical  parameters  of  the  transport  channel  and  the  sense  structure. 

Theoretical  Work 

Time  and  frequency  domain  responses  for  a  line  charge  traveling  in  a  single  tap  current  sens¬ 
ing  electrode  structure  (Fig.  5(a))  are  determined  by  computing  the  charge  induced  upon  the  sense 
electrodes  by  the  line  charge.*5*  Analytical  results  for  the  frequency  response  are  obtained  for  arbi¬ 
trary  charge  depth  when  the  gaps  between  electrodes  are  negligible,  while  numerical  results  are 
obtained  for  non-zero  gaps.  Normalized  time  and  frequency  response  data  for  the  structure  of 
Fig.  5(a)  for  several  values  of  charge  depth  and  gap  size  are  generated.  Figure  6  shows  the  time 
domain  response  for  the  structure  of  Fig.  5(a)  with  a  vanishingly  small  gap  between  the  surface 
olates.  The  spreading  of  the  response  about  the  center  is  greater  for  larger  charge  depth  due  to  the 
greater  spreading  of  lines  of  flux  emanating  from  the  charge  packet.  Hence,  the  rolloff  in  the  fre¬ 
quency  response  is  greater  for  larger  depths  of  transport.  The  presence  of  gaps  between  current 
sense  electrodes  introduces  additional  rolloff  in  the  frequency  response.  A  small  signal  equivalent 
circuit  for  current  sensing  consists  of  a  single  current  source  shunted  by  the  structure  capacitance. 

The  frequency  response  for  the  single  tap  voltage  sensing  structure  of  Fig.  5(b)  has  been  com¬ 
puted  analytically  in  terms  of  a  superposition  of  current  sense  responses.*5*  Voltage  sensing  intro¬ 
duces  a  multiplicative  factor  into  the  frequency  response,  which  is  dependent  upon  the  length  of 
the  sense  electrode.  It  is  found  that  a  null  in  the  frequency  response  may  be  ntroduced  by  proper 
choice  of  sense  electrode  length.  An  equivalent  circuit  for  the  single  tap  voltage  sensing  ructure 
consists  of  a  voltage  source  in  series  with  the  structural  capacitance  of  adjacent  electrodes.  It  is 
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Figure  5(a).  AC  model  of  simple  current  sensing  structure. 


Metal  Plate  A 

Metal  Plate  B 

ri 

y 

\ 

Li 

/ - ►)  |  (■* - 

Line  Charge  QL  9  1 


Figure  5(b).  AC  model  of  single  electrode  voltage  sensing  structure. 
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found  that  frequency  response  rolloff  is  lessened  for  shorter  electrode  length. 

The  effect  of  the  distribution  of  charge  in  the  packet  upon  current  and  voltage  sensing  fre¬ 
quency  responses  is  determined  numerically.'5'  using  previous  results  for  the  packet  charge  distri¬ 
bution  '4'  For  small  charge  loads,  the  rolloff  introduced  by  the  packet  shape  at  low  frequencies  is 
much  smaller  than  the  more  fundamental  current  and  voltage  sensing  frequency  response  rolloffs. 

Theoretical  work  has  also  focused  upon  nonlinearities  introduced  into  the  frequency  response 
of  the  ACT  device  due  to  the  finite,  non-point  nature  of  the  surface  acoustic  wave  (SAW)  as  a 
sampler  of  charge.  Figure  7  illustrates  the  dependence  of  the  normalized  charge  density  in  a  single 
charge  packet  as  a  function  of  time  for  50%  normalized  charge  load  as  it  flows  in  the  ACT  chan¬ 
nel.'5'  It  is  found  that  the  size  and  shape  of  the  packet  is  a  function  of  the  amount  of  charge  con¬ 
tained  in  the  packet.  Therefore,  the  spectrum  (Fourier  transfoim)  of  the  charge  packet  is  also  a 
function  of  charge  load.  This  gives  rise  to  the  introduction  of  new  frequencies  (harmonics)  in  the 
frequency  response  of  the  ACT  device  to  sinusoidal  steady  state  input. 

For  a  single  tone  input,  a  theoretical  method  has  been  devised  for  the  computation  of  the  level 
of  the  harmonic  components  present  in  the  output  signal  as  a  function  of  the  magnitude  and  fre¬ 
quency  of  the  input  tone.  Initial  computations  with  this  method  were  limited  to  cases  for  which 
the  signal  frequency  is  an  integer  fraction  of  the  SAW  frequency,  but  recently  a  method  has  been 
found  to  treat  the  general  case.  Shown  in  Fig.  8  is  a  plot  of  the  magnitude  of  the  second  harmonic 
versus  the  magnitude  of  the  fundamental  signal,  for  fSignai/fsaw  =  1/7  and  0  to  50%  signal  charge 
load  with  a  50%  quiescent  charge  load.  The  curve  shows  the  expected  quadratic  behavior,  since  this 
is  a  second  order  effect.  For  a  signal  charge  load  of  40%,  the  curve  indicates  that  second  harmonic 
supression  of  55dB  below  the  fundamental  component.  Computation  of  the  third  harmonic  for  this 


case  has  also  been  carried  out. 
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Future  Work 


Computation  of  harmonic  components  for  cases  with  arbitrary  relation  between  SAW  and  sig¬ 
nal  frequency  will  be  carried  out.  as  will  simulation  of  the  two-tone  input  signal.  Experiments 
will  be  designed  to  verify  the  theory.  Additionally,  theory  of  large  signal  detection  in  NDS  struc¬ 
tures  will  be  treated. 

Experimental  Work 

The  voltage  sensing  frequency  response  of  a  traveling  wave  charge  transport  device  whose  AC 
model  is  shown  in  Fig.  5(b)  is  measured  and  compared  with  theoretical  results.  The  experimental 
and  theoretical  frequency  responses  are  shown  in  Fig.  9.  The  agreement  between  theory  and  exper¬ 
iment  is  very  good  in  the  low  and  mid-frequency  range.  The  major  source  of  error  at  high  frequen¬ 
cies  is  due  to  uncertainity  in  charge  transport  depth.  It  is  calculated  that  this  uncertainty  can 
account  for  the  difference  between  theory  and  experiment. 

Conclusion 

An  analytic  and  numerical  investigation  of  the  frequency  responses  of  non-destructive  sensing 
structures  in  a  traveling  wave  charge  transfer  device  is  performed  and  is  shown  to  be  in  agreement 
with  experimental  data  taken  on  a  single  tap  voltage  sensing  structure.  The  theory  provides  a 
foundation  for  accurate  electrical  modeling  of  arbitrary  non-destructive  sense  structures. 
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3.  Transport  Research 


Objective 

The  objective  of  this  work  is  to  develop  the  theory  and  computer  programs  to  calculate  the 
two-dimensional  potentials  in  the  transport  channel  of  an  acoustic  charge  transport  (ACT)  device 
for  various  device  architectures  and  to  ascertain  the  effects  of  these  potentials  on  the  basic  tran¬ 
sport  properties  of  an  ACT  device. 

Overview 

There  are  several  tasks  of  significance  in  this  research.  One  task  is  to  calculate  the  electros¬ 
tatic  potentials  in  the  channel.  Another  task  is  to  analyze  the  transfer  inefficiencies  in  an  ACT  dev¬ 
ice  caused  by  surface  structures  such  as  gaps  between  Schottky  plates  and  voltages  on  electrodes  in 
the  channel.  Also  of  interest  is  charge  storage  in  an  ACT  device  which  corresponds  to  the  case  of 
applying  a  large  negative  voltage  to  the  electrodes  in  the  channel.  Finally,  there  is  the  t.i^k  of 
analyzing  charge  extraction  in  an  ACT  device  which  corresponds  to  the  extreme  case  of  applying  a 
large  positive  voltage  to  the  electrodes  in  the  channel. 

Potentials 

Figure  10  shows  typical  structures  for  which  the  potentials  are  desired.  Typically,  these 
structures  consist  of  an  n-type  active  layer  on  top  of  a  p-type  substrate,  or  a  semi-insulating  sub¬ 
strate  that  is  modeled  as  p-type  due  to  trapping  effects,  with  Schottky  plates  on  the  surface  of  the 
n-type  active  layer.  The  size  of  the  gaps  are  not  small  enough  to  be  neglected  and  must  be  included 
in  the  analysis. 

The  key  to  solving  for  the  desired  potentials  in  a  simple  fashion  is  to  separate  the  problem 
into  two  parts  and  to  utilize  superposition  to  obtain  the  total  solution.  Figure  11  illustrates  the 
separation  of  the  problem.  All  of  the  potential  variation  due  to  space  charge  is  separated  into  a 
one-dimensional  problem.  The  solution  to  this  problem  can  be  found  analytically  and  is  added  to 
the  solution  of  a  two-dimensional  Laplace  problem.  This  two-dimensional  problem  is  picked  such 
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Figure  10. 


Typical  structures  encountered  in  buried  channel  charge  transfer  devices. 

(a)  Isolated  gap  with  symmetrical  biasing. 

(b)  Isolated  gap  with  unsymmetricai  biasing. 

(c)  Isolated  sense  electrode. 

(d)  Arbitrary  periodic  structure. 
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Figure  11.  Separation  of  a  potential  problem  into  two  problems. 

(a)  Two-dimensional  potential  problem  with  space  charge. 

(b)  One-dimensional  potential  problem  containing  all  space  charge. 

(c)  Two-dimensional  potential  problem  without  space  charge. 
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that  subsequent  superposition  with  the  one-dimensional  potential  satisfies  all  differential  equations 
and  boundary  conditions.  This  two-dimensional  Laplace  problem  is  then  solved  using  Fourier 
analysis  or  conformal  mapping.  Full  details  of  these  calculations  are  in  reference  [6]. 

By  using  these  techniques  one  can  obtain  the  empty  channel  potentials  for  any  arbitrary 
structure  and  doping  profile.  Figure  12  shows  the  normalized  potential  deviation  from  the  regular 
one-dimensional  potential  caused  by  a  single  gap  with  gap  width  to  active  layer  ratio  of  2:3.  To 
obtain  the  actual  potential  deviation  for  a  given  doping  profile  and  channel  minimum  point  one 
must  multiply  the  normalized  potential  deviation  by  the  scaling  factor. 

where  d  is  the  active  layer  depth,  ym  is  the  transport  depth,  p( y)  is  the  charge  density  of  the  chan¬ 
nel.  is  the  permittivity  of  the  substrate,  and  e2  is  the  permittivity  of  the  dielectric  above  the 
substrate.  Plots  for  other  gap  dimensions  are  in  reference  [l]. 

The  above  analysis  doesn't  include  mobile  charge  in  the  channel  however,  it  is  still  useful  in 
analyzing  charge  storage,  charge  extraction,  and  transfer  inefficiency.  A  way  has  been  theoretically 
derived  to  include  mobile  charge  packets  in  the  channel  and  to  also  find  the  equilibrium  boundaries 
associated  with  these  charge  packs.  This  method  avoids  the  necessity  of  evaluating  the  two- 
dimensional  convolution  integral  previously  used  in  evaluating  charge  packet  potentials  and  also 
avoids  a  finite  difference  approach.  It  has  the  promise  of  greatly  reducing  the  computation  time 
required  to  model  mobile  charge  in  the  channel.  At  present  the  programming  required  is  not  fully 
completed  and.  hence,  no  results  are  currently  available. 

Transfer  Inefficiency 

The  desired  quantities  of  interest  in  this  area  are  the  transfer  inefficiencies  as  a  function  of  gap 
size,  transport  depth.  SAW  power  and  induced  or  applied  electrode  voltage.  The  first  step  in  calcu¬ 
lating  the  transfer  efficiency  is  to  calculate  the  reduction  in  charge  capacity  as  a  function  of  various 
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parameters.  A  one-dimensional  charge  packet  model  predicts  that  the  charge  capacity  reduction 
factor  (ccrf)  as  a  function  of  Ex,  which  is  the  x-component  of  the  perturbing  field,  is 


ccrf  = 


where  k  is  the  wave  vector  of  the  SAW. 

The  effect  of  gaps  in  the  channel  on  charge  transport  is  a  problem  of  great  concern.  Figure  13 
shows  the  normalized  peak  perturbing.  Ex,  as  a  function  of  normalized  gap  size  for  various  normal¬ 
ized  transport  depths.  The  parameter  g  is  the  gap  size  and  the  parameter  d  is  the  total  depletion 

depth  of  the  device.  To  unnormalize  the  result  one  multiplies  by  the  scaling  factor 
y  m 

SF2  =  -  J  p(y)dy.  The  ccrf  is  then  easily  found  for  various  SAW  frequencies  and  poten- 

el  +  €2  O 

tials.  If  the  maximum  operating  charge  capacity  of  an  ACT  device  exceeds  the  adjusted  charge 
capacity  given  by  the  maximum  SAW  capacity  less  the  percentage  given  by  the  ccrf,  then 
significant  transfer  inefficiency  can  be  expected;  otherwise,  the  transport  will  be  minimally  per¬ 
turbed.  A  similar  analysis  can  be  done  to  find  the  ccrf  due  to  small  voltages  on  electrodes  in  the 
channel. 


Charge  Storage 

In  analyzing  charge  storage  the  desired  quantities  of  interest  are  the  threshold  voltage  of 
storage  and  the  cutoff  voltage  of  transport.  Given  a  few  assumptions  and  approximations  it  is  pos¬ 
sible  to  generate  estimates  of  these  quantities. 

The  first  assumption  is  that  a  pulse  will  be  used  to  store  the  charge.  The  ramification  of  this 
is  that  the  deep  traps  in  the  semi-insulating  layer  will  not  react  with  the  barrier  pulse  due  to  the 
long  time  constants  of  the  traps.  Also,  if  the  length  of  the  barrier  electrode  is  long  enough  such 
that  it  can  be  assumed  to  be  semi-mfinite,  then  the  fields  due  to  two  coplanar.  semi-infinite,  surface 
plates  separated  by  a  gap  are  sufficient  to  describe  the  barrier  fields.  The  maximum  Ex  of  the  struc¬ 
ture  for  a  given  transport  depth  provides  the  stopping  force  to  charge  transport.  A  good  estimate 
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Normalized  Peak  Perturbing  Field 


Figure  13.  Perturbing  E-field  vs  gap  size. 
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of  the  maximum  Ex  is  Ex  down  the  center  of  the  gap.  For  large  transport  depths,  this  is  the  max¬ 
imum  Ex.  Since  the  fields  in  the  described  problem  scale  linearly  with  the  barrier  voltage  Vb  and 
the  gap  width  g,  a  normalized  plot  of  Ex  can  be  made  and  is  depicted  in  Fig.  14.  Fig  14  shows 
Enorm-  normalized  Ex,  versus  ynorm.  the  normalized  ymi,i.  The  following  relations  are  used  to  find  Ex 
for  a  given  ymill: 

Ex  —  Exnorm  Vb/g 

y  min  y  norm? 

The  threshold  voltage  of  storage  is  achieved  when  Ex  equals  the  product  of  Vsaw  and  k.  The  cutoff 
voltage  of  transport  is  somewhat  higher,  for  to  guarantee  that  no  charge  leaks  past  the  barrier  the 
Ex  of  the  barrier  must  overcome  the  peak  SAW  E-field,  the  induced  E-field  of  the  stored  charge,  Ein, 
and  the  barrier  must  also  supply  a  sufficient  potential  height  on  top  of  this  to  prevent  diffusion 
over  the  barrier.  Ein  can  be  approximated  by  an  infinite  train  of  line  charges  simulating  the  packets 
stored  behind  the  barrier.  Table  1  gives  estimates  of  maximum  Em  for  various  ratios  of  ymin  to 
SAW  wavelength.  Thus.  Vb  can  be  determined  for  the  total  required  Ex  from  Fig.  14.  In  addition, 
to  ensure  a  barrier  against  diffusion,  twice  the  desired  diffusion  barrier  should  be  added  to  the 
result  to  predict  the  cutoff  voltage  of  transport.  This  is  because  the  first  half  goes  to  raising  the 
potential  of  the  stored  charge  and  the  second  half  goes  towards  stopping  the  diffusion. 

Charge  Extraction 

This  phenomenon  is  by  far  the  most  difficult  to  analyze.  The  extract  voltage  and  pulse  shape 
of  the  extracted  packet  depend  on  the  doping  profile  as  well  as  the  structure.  To  extract  charge 
from  the  channel  it  is  necessary  to  have  a  clear  trajectory  to  the  extract  electrode.  In  general,  the 
shallower  transport  depths  require  less  extraction  voltage.  Table  2  yields  the  predicted  extract  vol¬ 
tages  with  reference  to  the  Schottky  plates  for  various  yinm  f°r  a  -  micron  Schottky  electrode 
between  2  micron  gaps  on  a  5  micron  active  layer  of  doping  101’  cm  \ 


Table  1.  Normalized  Induced  E-field 


Y  min 
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.1 
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Em 


.0996 

.338 

.613 

.865 

1.079 


£  _  Qmax 
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W  =  width  of  channel 
Qmax  ===  maximum  charge  capacity  of  a  well 
=  induced  E-field 


Ex(total)  =  VSAW  *  k  +  E^ 


Table  2.  Extract  Voltages 


ymin 

extract 

2.5 

9.395 

2.0 

8.456 

2.25 

7.516 

1.75 

6.577 

1.5 

5.637 

4.  Low  Process  Temperature  Dielectric  Research 


Objective 


The  objective  of  this  research  is  to  develop  a  low  temperature  fabrication  process  for  mul¬ 
tilevel  interconnections  using  a  polvimide  dielectric.  Material  properties  of  the  dielectric  processed 
with  the  above  low  temperature  fabrication  procedure  are  also  studied. 
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Experimental  Work 


Development  of  the  wet  chemical  fabrication  process  includes  polyimide  dielectric  selection, 
polyimide  etchant  evaluation,  and  comprehensive  parametric  optimization  of  basic  processing 
parameters.  Figure  15a  and  15b  outline  the  basic  steps  in  the  fabrication  process. 

Lnique  qualities  of  the  interlayer  dielectric  and  fabrication  process  are: 

a.  Low  dielectric  constant  to  reduce  capacitive  coupling  between  adjacent  metaliza- 
tions. 

b.  Application  of  two  polyimide  layers  to  achieve  low  pinhole  density,  reducing  inter¬ 
level  shorts. 

c.  Etchant  selection  which  achieves  minimal  degradation  of  the  GaAs  substrate  and 
first-level  metal  surface  morphology  during  processing. 

d.  Low  temperature  polyimide  cure  to  eliminate  alloying  of  the  first-level  Schottky 
barrier  features  into  the  GaAs  substrate. 

e.  Favorable  edge  profile  for  uninterrupted  metalizations  while  providing  small 
geometry  capability.  See  Figs.  15c  through  15f. 


The  electrical  and  mechanical  properties  of  polyimide  films  are  evaluated  using 
three  test  vehicles:  capacitor,  via,  and  crossover  arrays.  The  following  polyimide 
film  characteristics  are  evaluated  and  corresponding  values  listed  in  Table  3:  dissi¬ 
pation  factor,  dielectric  constant,  dielectric  breakdown  strength,  pinhole  density, 
and  via  resistance. 

The  completeness  ot  the  low  temperature  polyimide  cure  can  be  determined  from  dissipation 
factor  measurements.  The  average  dissipation  factor  of  capacitors  with  plate  areas  (809  /xm)2,  (572 
/xm)“.  and  (404  ,um)2  is  similar.  This  similarity  of  dissipation  factors  suggests  a  bulk  limited  loss 
mechanism. 


The  relative  dielectric  constant  is  calculated  from  capacitance  data  taken  during  dissipation 
factor  measurements. 

Dielectric  breakdown  strength  is  measured  using  a  capacitor  test  structure.  Voltage  across  tiie 
capacitor  plates  is  slowly  increased  until  the  dielectric  fails. 
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First  Polyimide 
Application 


Second  Polyimide 
Application 


Photoresist 

Layer 


Contact  Mask 


Process  sequence  for  electrical  connection  between  points  A  and  C  without  connection 
to  point  B. 

a.  Substrate  with  first-level  metalization  present. 

b.  Apply  first  polyimide  coating  and  soft-bake. 

c.  Apply  second  polyimide  coating  and  soft-bake, 
d  Apply  photoresist  and  soft-bake. 

e.  Align  and  expose. 
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Table  3 


dissipation  factor  0.0143  ±  0.0004  «(  1  MHz 

Dielectric  constant  between  3.9  and  4.2  at  1  MHz 

Dielectric  breakdown  strength  193  ±  60  V/u  m  DC 
Pinhole  density  less  than  1.17/cm2 

V  a  resistance  between  2  and  6  ohms 

The  total  area  tested  for  interlevel  shorts  from  pinholes  measured  85  mm2.  Crossover  arravs 
are  used  for  this  measurement.  No  shorts  were  found. 

Resistance  of  vias  with  contact  areas  of  (100  fim)2.  (50  yam)2,  and  (10  /am)2  is  measured.  No 
correlation  between  via  contact  area  and  via  resistance  is  noted.  It  is  believed  that  only  point  con¬ 
tact  between  adjacent  metalizations  is  achieved  due  to  the  oxide  of  the  first-level  metal.  Lower 
resistance  could  probably  be  obtained  if  the  oxide  of  the  first-level  metal  could  be  more  effectively 
removed,  or  if  vias  are  sintered.  However,  sintering  would  alloy  Schottky  barrier  features  into  the 
GaAs  substrate. 

A  more  thorough  discussion  of  the  dielectric  processing  can  bo  found  in  reference  [7], 
Conclusion 

A  low  temperature  process  to  "attern  polyimide  dielectric  is  achieved.  Polyimide  has  ade¬ 
quate  electrical  and  mechanical  properties  for  interlayer  isolation.  The  fabrication  procedure 
developed  is  compatible  with  both  photolithographic  and  process  parameter  constraints.  Improve¬ 
ments  in  reducing  via  contact  resistance  need  to  be  addressed. 
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5.  Deep  Implantation  Research 


Objective 

The  overall  objective  of  this  research  is  to  study  the  fundamental  properties  of  deep  implants 
into  GaAs  and  explore  the  basic  processing  techniques  necessary  to  apply  deep  implants  to  GaAs 
ICs.  Work  to  date  has  concentrated  on  one  specific  form  of  deep  implant  into  GaAs:  the  isolation 
implant  The  purpose  of  this  form  of  implant  is  to  create  a  disordered  layer  in  the  crystal  struc¬ 
ture,  thereby  greatly  increasing  the  resistivity  of  the  material.  This  technique  is  very  useful  for 
providing  electrical  isolation  between  active  regions  on  a  GaAs  substrate. 

Studies  thus  far  have  been  in  two  areas.  The  first  of  t,  ase  is  the  development  of  GaAs  compa¬ 
tible  processing  techniques.  In  this  area,  the  primary  goal  has  been  the  development  of  a  photol- 
ithographically  definable  deep  implant  mask  capable  of  completely  stopping  the  light  high  energy 
particles  used  in  isolation  implants.  The  second  area  of  study  has  been  research  into  the  fundamen¬ 
tal  aspects  of  deep  implants.  Here,  work  is  being  done  on  examining  the  general  characteristics  of 
isolation  implants  both  into  bare  GaAs  and  through  thin  films  on  GaAs.  A  comparison  of  these 
general  characteristics  will  hopefully  lead  to  more  effective  masking  strategies  and  better  device 
isolation.  Since  studies  in  this  area  will  be  done  using  both  proton  and  deuteron  implant  isolation, 
comparisons  of  the  basic  properties  of  each  of  these  will  also  be  obtained. 

Deep  Implant  Mask 

Several  criterion  exist  which  the  deep  implant  mask  must  meet  in  order  to  fulfill  both  the 
requirement  of  complete  masking  and  the  needs  of  specific  applications.  The  mask  must  be  photol- 
ithograp’  ically  definable,  have  high  resolution  and  high  aspect  ratios,  be  thick  enough  to  com¬ 
pletely  stop  light,  high  energy  particles,  and  oe  easily  removable.  These  criterion  can  not  be  met 
using  conventional  implant  masking  processes  due  to  the  low  mass  and  high  energy  of  the  incident 
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One  masking  technique  capable  of  achieving  these  goals  which  appears  particularly  promising 
has  been  the  focus  of  this  investigation.  A  cross  sectional  view  of  this  mask  appears  in  Hg.  16 
The  purpose  of  each  of  the  layers  is  as  follows: 

***  Dielectric  Layer:  The  primary  function  of  this  layer  is  to  allow  easy  removal  of  the  mask 
after  the  high  energy  implant  has  been  performed. 

***  Thin  Metal:  This  layer  functions  as  an  electrode  during  the  plating  of  the  thick  metal  layer. 

***  Thick  Photoresist:  This  layer  confines  the  plating  to  specific  areas  and  prevents  lateral 
spreading  of  the  plated  metal.  It  is  removed  prior  to  implantation. 

«**  Thick  Metal:  This  metal  stops  incident  high  energy  particles,  masking  the  substrate  in  areas 
where  it  is  plated. 

This  mask  has  been  effectively  used  for  protecting  areas  of  thick  epitaxial  layers  on  GaAs 
during  proton  isolation.  Preliminary  results  are  very  promising,  indicating  that  masked  areas  were 
well  protected  from  the  implant. 

Isolation  Implant  Research 

Although  this  work  is  in  the  early  stages,  some  preliminary  investigations  have  been  per¬ 
formed.  Most  noteworthy  of  these  has  been  the  investigation  of  epitaxial  layer  resistivity  change 
with  implant  isolation  dose  for  295  KeV  protons.  Using  a  substrate  prepared  as  in  Fig.  17,  the 
resistance  of  the  epitaxial  layer  was  measured  as  the  proton  dose  was  incrementally  increased. 
Using  measured  length  and  width  values  for  the  substrate,  these  resistance  measurements  were 
then  used  to  calculate  the  sheet  resistance  of  the  epitaxial  material.  The  resultant  graph  of  sheet 
resistance  versus  proton  dose  appears  in  Fig.  18.  The  inset  graph  shows  the  doping  profile  of  the 
substrate  used  in  this  experiment.  From  this  graph,  some  general  conclusions  about  the  properties 
of  proton  isolation  implants  may  be  drawn.  First,  the  shape  of  the  graph  indicates  that  carrier 
removal  per  incident  ion  drops  near  zero  above  a  certain  value  of  sheet  resistance.  In  these  experi- 
menis.  this  value  was  found  to  be  approximately  400  megaohms  per  square.  Assuming  an  epitaxial 
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layer  thickness  of  2.75  microns,  this  corresponds  to  a  material  resistivity  of  110  kilo-ohm-cm. 
This  value  is  somewhat  lower  than  values  previously  reported  in  the  literature.  Two  effects  may 
be  responsible  for  this.  The  first  effect  is  bulk  conduction  through  the  substrate.  This  effect 
becomes  important  as  the  material  resistivity  of  the  proton  bombarded  epitaxial  layer  approaches 
the  resistivity  of  the  semi-insulating  substrate,  the  second  effect,  surface  leakage,  also  produces  a 
shunt  path  across  the  epitaxial  layer.  Once  again,  this  will  place  an  upper  limit  on  measurable 
material  resistivi  y. 

Another  useful  result  which  may  be  extracted  from  this  data  set  is  the  average  carrier  remo¬ 
val  per  incident  proton.  This  number  is  not  intended  to  serve  as  a  tool  whereby  exact  dosages  tc 
isolate  epitaxial  layers  can  be  calculated,  since  carrier  removal  has  been  found  by  previous 
researchers  to  varv  with  depth.  Rather,  it  allows  one  to  get  a  macroscopic  view  of  the  carrier 
removal  process.  To  do  the  calculation,  sheet  resistance  is  first  converted  into  ns.  the  average  free 
carrier  concentration  per  cm",  using  the  relationship  ns  =  (q/cR.)  T  As  increasing  doses  are  applied 
to  the  substrate,  ns  decreases.  The  change  in  free  carrier  concentration  divided  by  the  change  in 
dose,  Ans/Adose.  yields  the  approximate  average  carrier  removal  per  proton.  In  the  calculation, 
mobi'ity  is  assumed  to  be  constant.  This  assumption  should  be  fairly  valid  since  mobility  is  deter¬ 
mined  by  impurity  concentration,  and  this  does  not  change  significantly  for  the  proton  doses  used 
in  the  experiment.  If  mobility  does  change  to  a  non-trivial  degree,  the  numbers  calculated  for  car¬ 
rier  removal  will  be  somewhat  higher.  For  the  data  set  presented  in  Fig.  18,  the  initial  value  of 
carrier  removal  at  low  doses  was  found  to  be  4.3  carriers/proton.  At  higher  doses,  this  value 
declined  significantly.  This  value  is  very  likely  to  be  low  since  the  range  of  protons  in  GaAs  is 
greater  than  the  thickness  of  the  epitaxial  layer  used  in  this  experiment.  Hence,  not  all  of  the  pro¬ 
tons  take  place  in  the  carrier  removal  process.  In  a  similar  experiment  performed  on  epitaxial 
lavers  having  a  thickness  greater  than  the  range  ol  protons  in  GaAs,  average  carrier  removal  rate 
was  found  to  vary  from  10.4  at  low  doses  to  0.14  carriers/proton  at  higher  doses.  This  data  is 
presented  in  Table  4. 
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Table  4 

Carrier  Removal  per  Proton  for  295  KeV  Isolation  Implant 
of  a  Thick  Epitaxial  Layer  on  GaAs 


Sheet  Resistance 
(ohms/square) 

ns  (cm  2) 

Proton  Dose 
cm-2 

Ans/Adose 

2.75  k 

3.79  E  11 

0 

0 

3.21  k 

3.25  E  11 

5.2  E  9 

10.38 

3.66  k 

2.85  E  11 

10.4  E  9 

7.69 

4.28  k 

2.43  E  11 

15.6  E  9 

8.08 

5.49  k 

1.90  E  11 

2.60  E  10 

5.1 

8.31  k 

1.25  E  11 

7.60  E  10 

1.3 

9.38  k 

1.11  E  11 

12.8  E  10 

0.27 

10.69  k 

9  74  E  10 

23.0  E  10 

0.14 

A  final  conclusion  which  may  be  drawn  from  the  evidence  to  date  is  that  isolation  also  occurs 
near  the  surface  of  the  epitaxial  layer.  This  conclusion  is  based  on  the  observation  that  fairly  high 
isolation  sheet  resistances  are  obtained  in  spite  of  doping  near  the  surface.  If  isolation  were  not 
occuring  near  the  surface,  much  lower  sheet  resistances  would  be  anticipated.  Because  isolation  is 
theorized  to  result  from  protons  inducing  disorder  in  the  crystal  lattice,  and  because  this  disorder  is 
usually  induced  near  the  end  of  an  individual  proton's  range,  isolation  near  the  surface  appears  at 
first  to  be  an  anomaly  if  one  assumes  a  simple  gaussian  model  for  the  distribution  of  protons  on 
GaAs.  The  work  to  be  described  next  offers  one  explanation  of  this. 
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Future  work  in  this  area  will  study  deuterium  isolation  and  will  involve  in  situ  monitoring 
of  the  change  in  doping  profile.  This  work  should  provide  sufficient  data  to  more  accurately  calcu¬ 
late  the  carrier  removal  per  ion  and  analyze  the  carrier  removal  versus  depth.  Previous  research1 151 
indicates  that  resistivity  for  deuterium  may  increase  to  a  plateau  and  then  actually  decrease  some¬ 
what  as  the  mechanism  of  hopping  conduction  begins  to  manifest  itself. 

A  second  investigation  which  has  yielded  valuable  data  studies  the  physical  distribution  of 
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implanted  protons  and  deuterons  in  GaAs.  In  this  study,  bare  GaAs  is  bombarded  with  one  of  the 
above  ions  at  an  energy  of  100.  200,  or  295  KeV  and  a  dose  of  2  X  1015  ions/cm2.  Next,  using 
secondary  ion  mass  spectrometry  (SIMS),  the  density  distribution  of  the  implanted  species  is 
observed.  This  is  easily  done  for  deuterium,  but  cannot  be  as  easily  accomplished  for  hydrogen  due 
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Figure  19. 


Distribution  of  295  KeV  deuterium  in  bare  [100]  cut  GaAs. 
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Table  5 


Range  of  Protons  and  Deuterons  in  Bare  GaAs 
and  Through  Thin  Surface  Films  on  GaAs 


100 

KeV 

200 

KeV 

295 

KeV 

Proton  Range  in  Bare  GaAs 

— 

— 

2.51 

Deuteron  Range  in  Bare  GaAs 

0.99 

1.75 

2.50 

Deuteron  Range  in  GaAs  Through 
1/2  micron  of  Polyimide 

— 

— 

2.0 

Deuteron  Range  in  GaAs  Through 

1  micron  of  Polyimide 

_ 

_ _ 

1.45 

*  All  Values  in  table  are  in  microns  * 

An  approximate  value  for  average  carrier  removal  per  proton  at  low  doses  was  found  to  be 
7-10  carriers/proton  for  a  thick  epitaxial  layer  doped  at  1013  cm-3 

Initial  studies  of  proton  and  deuteron  distribution  in  GaAs  as  a  result  of  ion  implantation  at 
energies  of  100,  200.  and  295  KeV  have  been  performed.  These  studies  have  suggested  an  explana¬ 
tion  for  the  observed  surface  region  isolation  in  spite  of  large  ion  range. 

Studies  of  deuterium  bombardment  through  surface  films  have  demonstrated  that  a  295  KeV 
deuteron  implant  is  attenuated  by  approximately  120  KeV  per  micron  of  PI-2555  polyimide. 


6.  Schottky  Barrier  Leakage  Research 


Objective 

The  objective  of  this  research  is  to  investigate  the  factors  which  influence  reverse  bias  leakage 
current  in  GaAs  Schottky  barrier  junctions.  This  work  is  accomplished  by  studying  the  effect  of 
diode  edge  guarding,  barrier  metal  type,  and  process  temperature  on  leakage  current  behavior 

Experimental  Work 

A  Schottky  barrier  diode  array  photomask  is  designed  to  permit  the  fabrication  of  diodes  with 
several  sizes  of  active  area.  The  mask  is  designed  so  that  half  of  the  diodes  are  constructed  with  a 
polyimide  guard  ring  to  reduce  edge  leakage,  while  the  other  half  are  unguarded. 

Diode  arrays  are  fabricated  with  different  barrier  metal  types  using  standard  GaAs  compati¬ 
ble  lift  off  processes  for  both  the  ohmic  contact  and  Schottky  level  pattern  steps.  The  photoresist 
developer  used  in  this  process  automatically  strips  the  surface  oxide  in  the  diode  active  area  before 
metallization  occurs  so  that  the  undesirable  effects  of  surface  oxide  are  minimized.  After  all  pho¬ 
tolithography  is  completed,  the  aiode  arrays  are  diced  into  small  groups,  subjected  to  well  con¬ 
trolled  heat  treatment  (annealing)  using  a  nitrogen  ambient  in  the  range  of  25°C  -  275°C,  and  DC 
probed  to  obtain  Schottky  barrier  height  and  reverse  leakage  current  data. 

Experimental  Results 

The  effect  of  polyimide  guard  rings  which  reduce  edge  leakage  is  shown  in  Fig.  20  where  the 
reverse  leakage  current  as  a  function  of  diode  area  is  plotted  for  constant  reverse  bias  voltage.  The 
theoretical  leakage  current  obtained  using  thermionic  emission  theory  with  the  measured  Schottky 
barrier  height  is  also  shown  for  comparison.  This  data  shows  the  relative  contributions  of  bulk 
and  edge  leakage  encountered  in  typical  GaAs  Schottky  barrier  diodes  of  this  size  since  the  leakage 
current  difference  between  the  guarded  and  unguarded  diodes  may  be  presumed  to  be  due  to  edge 
leakage.  Although  the  barrier  in  this  experiment  is  a  multi-layer  combination  of  Cr  ( 1 00 
Angstroms)  and  Cu/Al  (1500  Angstroms),  the  Cr  metal  determines  the  diode  characteristics 
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because  it  is  the  metal  layer  which  contacts  the  GaAs  surface.  The  edge  leakage  component  derived 
from  tests  like  these  can  vary  significantly  from  one  batch  of  diodes  to  the  next  and  in  some  cases 
it  can  be  an  order  of  magnitude  higher  than  the  corresponding  bulk  value.  This  variation  is 
believed  to  be  due  to  varying  surface  contamination  levels  which  are  very  process  dependent. 
Although  variations  in  measured  Schottky  barrier  heights  are  obtained,  once  the  barrier  height  is 
experimentally  determined,  the  theory  predicts  the  experimental  leakage  very  accurately  as  shown 
in  Fig.  20.  In  general,  the  guarded  diode  leakage  is  always  slightly  greater  than  the  theoretical 
prediction.  Experiments  indicate  that  this  is  due  to  the  presense  of  generation-recombination  leak¬ 
age  components  which  are  not  accounted  for  in  the  barrier  leakage  theory. 

The  impart  of  barrier  metal  type  on  leakage  current  is  shown  in  Fig.  21.  where  leakage 
current  is  plotter  as  a  function  of  diode  area  for  unguarded,  unannealed  diodes  wit  v  barrier  metal 
type  as  a  parameter.  Barrier  height  measurements  show  that  the  metals  with  the  lowest  leakage 
currents  also  have  the  largest  barrier  height.  Hence,  the  var.ation  in  unannealed  leakage  current 
with  metal  type  is  directly  attributable  to  the  varying  work  functions  (and  hence  barrier  heights) 
which  the  metals  possess. 

Schottky  diodes  which  are  subjected  to  heat  treatment  can  exhibit  large  changes  in  1-akage 
current.  An  example  of  this  phenomenon  is  shown  in  Fig.  22  where  the  leakage  current  in  A1 
diodes  is  plotted  as  a  function  of  reverse  bias  with  annealing  temperature  as  a  parameter,  for  an 
anneal  lime  of  60  minutes.  A  leakage  current  reduction  of  roughly  two  orders  of  magnitude  is 
obtained  b”  annealing  the  A1  diodes  at  275°C  for  one  hour.  This  effect  has  been  observed  by  other 
workers*3'  and  is  due  to  the  formation  of  a  thin  AlGaAs  layer  at  the  metal  semiconductor  interface 
which  raises  the  effective  barrier  height  of  the  junction.  The  precise  effect  of  annealing  is  dependent 
on  the  particular  metal  and  the  temperature.  Both  leakage  current  increases  and  reductions  have 
been  observed  for  various  metals. 
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Figure  21. 


Leakage  current  versus  diode  area 


for  three  different  Schottky  metallizations 
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7a.  GaAs  SAW  Waveguide  Research 
Objective 

The  objective  of  this  work  is  to  study  SAW  generation  and  propagation  in  GaAs  acoustic 
waveguides.  Part  of  the  work  involves  gaining  more  detailed  knowledge  of  losses  in  basic  SAW 
propagation  on  GaAs,  and  energy  loss  mechanisms  in  the  surface  wave  transducers.  The  remainder 
involves  investigation  of  the  use  of  beam  compressors  and  waveguides  to  focus  and  confine  modest 
amounts  of  total  SAW  power  to  achieve  large  amplitudes  in  narrow  channel  acoustic  systems.  One 
embodiment  to  be  studied  involves  an  acoustic  resonant  ring:  a  waveguide  "racetrack"  which  reuses 
the  acoustic  energy  for  as  long  as  possible. 

Theoretical  Work 

A  mathematical  model  has  been  derived  in  order  to  study  losses  in  SAW  transducers  having 
internal  reflections,  Loss  terms  dealing  with  finger  resistance  and  SAW  propagation  losses  have 
been  derived  and  added  to  the  differential  equations  of  an  existing  coupled-mode  model  by  Hart¬ 
mann  et  alJ9*  The  solutions  to  the  new  differential  equations  have  been  used  in  a  computer  model¬ 
ing  program. 

The  program  calculates  SAW  amplitude  and  electrical  drive  current  as  a  function  of  either 
position  in  the  transducer  or  frequency.  Outputs  from  this  program  agree  well  with  experimental 
data  on  reflective  transducers  on  ST  quartz,  YZ  LiNb03,  and  GaAs. 

Our  model  predicts  that  finger  resistance  losses  are  more  serious  than  propagation  loss,  even 
though  the  internal  reflections  cause  the  SAW  energy  to  bounce  around  inside  the  transducer 
longer.  This  would  imply  that  the  metal  thickness  of  transducer  fingers  should  be  increased. 

Other  theoretical  work  involving  waveguide  calculations  such  as  mode  shapes  and  velocities 
uncovers  an  important  detail.  The  pseudo-SAW  has  propagation  loss  which  depends  strongly  on 
the  direction  in  which  the  wave  is  launche  l0*  (Fig.  23)  The  loss  is  minimum  exactly  aiong  a 
[llOj  direction,  and  increases  sharply  at  small  angles  up  to  about  6  degrees  awav.  A  waveguide 
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mode  is  equivalent  to  a  set  of  rays  which  zigzag  along  the  guide  at  some  angle  to  the  guide  axis,  and 
this  angle  decreases  with  increasing  guide  width.  Thus,  propagation  loss  should  vary  inversely 
with  guide  width.  We  have  recently  verified  this  experimentally. 

Experimental  Work 

A  high-precision  automated  laser  probe  has  been  built  for  this  work.  Its  demonstrated  capa¬ 
bilities  include  a  frequency  range  greater  than  1  GHz.  position  resolution  of  0.1  /xm.  minimum  laser 
spot  size  of  2  fj. m.  large  dynamic  range,  .id  amplitude  repeatability  of  about  0.3%.  A  block 
diagram  of  the  laser  probe  system  is  shown  in  Fig.  24. 

The  probe  has  been  used  to  measure  SAW  propagation  data  in  semi-insulating  GaAs.  which  is 
summarized  in  Table  6.  It  should  be  noted  that  the  material  measured  so  far  is  not  exactly  (100) 
cut,  but  the  more  common,  material  whose  plate  normal  is  rotated  two  degrees  away  from  the  [100] 
direction.  Some  crystal  growers  need  this  for  growth  of  good  epitaxial  layers.  These  measurements 
will  be  repeated  for  several  different  crystal  orientations,  including  exact  { 100}  cut. 

Conclusion 

Several  important  observations  may  be  made  from  this  data: 

1)  Unguided  propagation  loss  rises  as  the  1.5  power  of  frequency,  as  observed  by  Melloch 
&  Wagers.'11'  We  cannot  explain  at  this  time  why  our  observed  losses  are  higher  than 
theirs,  but  feel  our  measurement  method  is  more  direct  and  hopefully  more  accurate 
than  theirs. 

2)  The  velocity  dispersion  under  an  aluminum  layer  is  stronger  than  predicted  by  Datta 
and  HunsingerJ12' 

3)  Guided  propagation  loss  does  increase  with  decreasing  guide  width,  as  predicted. 


58 


Table  6. 


Measurements  of  Guided  &  Unguided  SAW  Propagation  Constants 
on  { 100}  Cut,  [l  10]  Propagating  GaAs 

I.  Unguided  Propagation 


360 

720 

1070 

MHz 

MHz 

MHz 

A.  Propagation  loss, 
dB/usec 

Metallized  surface 

.971 

3.4 

5.28 

Free  surface 

.985 

— 

5.58 

B.  Velocity,  meters/sec  [*] 

Metallized  surface 

2873.1 

2885.9 

2904.3 

+/-0.1 

+/-0.1 

+/-4.9 

Free  surface 

2861.5 

2862.9 

2859.8 

+/-0.2 

+/-3.5 

+/-0.8 

II.  Guided  Propagation  (delta-V/V  slot  waveguide) 


360 

1070 

MHz 

MHz 

Propagation  Loss,  dB/usec 

===== 

=== 

Wide  Guide  (10  lambda  @  360) 

1.59 

7.38 

Narrow  Guide  (3  lambda  @  360) 

3.04 

8.61 

B.  Velocity,  meters/sec  [*] 


Wide  guide,  fundamental 

2871.4 

symmetric  mode 

+/-0.2 

Narrow  guide,  fundamental 

2884.9 

symmetric  mode 

+/-0.4 

Narrow  guide,  first 

2917.0 

anti-symmetric  mode 

+/-0.4 

[*]  Note:  The  velocity  tolerances  stated  above  are  preliminary,  and  are  only  meant  to  give  some 
idea  r  l  the  relative  uncertainties  of  the  measurements. 
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7b.  GaAs  SAW  Propagation  Research 


Objective 

The  objective  of  this  research  is  to  investigate  the  details  of  SAW  propagation  in  GaAs.  This 
includes  the  development  of  models  which  can  be  used  to  predict  the  propagation  of  beam  profiles 
through  bare  GaAs  and  through  metal  grating  structures.  The  work  will  also  include  the 
verification  of  these  models  by  laser  probe  experiments. 

Theoretical  Work 

A  computer  program  has  been  developed  which  uses  standard  diffraction  theory'12'  to  predict 
the  beam  profile  propagated  a  distance  specified  by  the  user.  This  program  allows  one  to  use  either 
an  experimentally  determined  slowness  surface  or  one  which  is  theoretically  derived.  To  illustrate, 
the  user  enters  into  the  program  a  laser  probe  scan  taken  at  position  (a)  shown  in  Fig.  25,  and  this 
scan  is  propagated  the  distance  specified  by  the  user  according  to  the  selected  slowness  surface.  A 
comparison  between  the  beam  profile  measured  with  the  laser  probe  at  position  (b)  and  the  profile 
predicted  for  position  (b)  by  the  diffraction  theory  program  is  presented  in  Fig.  26.  For  this  calcu¬ 
lation  an  isotropic  slowness  surface  was  assumed  and  it  is  clear  from  the  inaccuracies  in  the  fine 
structure  of  the  beam  profile  that  the  isotropic  slowness  surface  does  not  closely  approximate  the 
GaAs  slowness  surface.  Because  no  experimental  data  is  available  in  the  literature  on  the  slowness 
surface  of  GaAs  a  method  of  computing  it  from  laser  probe  scans  was  developed.  The  details  of 
this  method  are  described  in  a  recently  completed  publication.'13' 

Experimental  Work 

Laser  probe  measurements  were  made  of  a  beam  propagating  away  from  a  bidirectional  trans¬ 
ducer  on  a  [l()0]-cut  Hoiizontal  Bridgman  Cr-doped  GaAs  substrate  and  the  slowness  surface 
obtained  is  shown  in  Figure  27.  The  velocity  and  attenuation  data  obtained,  compares  favorably 
with  theoretical  predictions  made  by  Penunuri191  and  the  measured  data  of  Deacon  and  Heigh  way 
which  are  also  presented  in  Figure  27. 
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Figure  26(a).  Measured  beam  profile  at  position  (b) 
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Velocity  Curves 


Penunuri  and  Lakm 
Deacon  and  Heighway 
Zero -Degree  Wafer 


This  experimentally  derived  slowness  surface  has  been  used  in  the  diffraction  theory  program 
and  very  accurate  predictions  of  both  the  magnitude  and  phase  of  the  beam  profile  have  been 
obtained  as  can  be  seen  in  Figure  28.  The  preciseness  of  the  beam  profile  predictions  validates  the 
theoretical  basis  of  the  diffraction  program  as  well  as  the  accuracy  of  the  slowness  surface  so  calcu¬ 
lated.  In  addition  to  these  measurements,  laser  probe  scans  were  made  for  a  range  of  kF  drive 
powers  and  no  appreciable  change  in  beam  profile  shape  was  observed.  This  preliminary  work  indi¬ 
cates  that  the  diffraction  oi  a  beam  propagating  in  unmetallized  GaAs  does  not  vary  with  acoustic 
power. 

Conclusion 

A  technique  for  computing  the  slowness  surface  of  GaAs  has  been  developed  and  the  results 
compare  favorably  with  the  theoretical  predictions  of  Penunuri.191  The  true  test  of  the  slowness 
surface  measurement  is  in  its  ability  to  predict  the  details  of  beam  profile  propagation  and  the 
slowness  surface  presented  herein  yielded  very  good  agreement  between  the  predicted  profile  and 
the  actual  profile.  In  addition,  laser  probe  scans  have  been  made  at  a  succession  of  transducer  RF 
drive  powers  and  no  appreciable  change  in  bea  diffraction  is  observed. 
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8.  Research  Summary 


Device  Physics  and  Modeling 

A  physical  model  for  traveling  wave  charge  injection  and  a  numerical  technique  for  solving 
this  model  have  been  established.  Preliminary  results  are  given. 

A  comprehensive  small  signal  theory  describing  the  frequency  response  and  sensitivity  of 
non-destructive  charge  sensing  electrodes  has  been  developed.  The  agreement  with  experimental 
data  is  excellent.  Future  work  in  this  area  will  focus  on  a  theoretical  description  of  large  signal 
operation. 

A  general  technique  for  determining  the  impact  of  arbitrary  static  surface  boundary  condi¬ 
tions  on  the  channel  potential  in  an  acoustic  charge  transport  (ACT)  device  has  been  developed. 
Simple  one-dimensional  models  are  given  for  charge  storage  and  transfer  inefficiency  due  to  gaps. 
Charge  extraction  voltages  can  be  calculated  from  the  empty-channel  potentials  and  examples  are 
given. 

Fabrication  Technology  Research 

A  fabrication  process  has  been  developed  for  denning  polyimide  dielectric  features  using  low 
process  temperatures.  The  electrical  characteristics  of  the  resulting  dielectric  have  been  evaluated 
and  found  to  give  adequate  DC  and  RF  characteristics  for  GaAs  integrated  circuit  and  crossover 
applications. 

An  initial  investigation  of  the  impact  of  diode  edge  guarding,  barrier  metal  type  and  process 
temperature  on  the  reverse  biased  leakage  current  of  Schottky  barrier  GaAs  diodes  has  been  per¬ 
formed.  The  work  indicates  that  drastic  reductions  in  leakage  current  may  be  obtained  through 
proper  choice  of  barrier  metal  type  and  process  temperature. 

A  fabrication  process  for  achieving  a  photolithographicallv  definable  high  resolution  mask 
which  is  capable  of  withstanding  high  energy  particle  bombardment  has  been  developed.  This  basic 
process  has  served  as  the  foundation  for  constructing  traveling  wave  charge  transier  device 


architectures  based  on  high  energy  proton  isolation.  Initial  studies  of  device  isolation  effectiveness 
have  been  completed. 

GaAs  Surface  Acoustic  Wave  Research 

GaAs  SAW  propagation  losses  in  unguided  and  guided  structures  have  been  accurately  meas¬ 
ured.  The  operation  of  GaAs  SAW  waveguide  structure  has  been  demonstrated.  A  methodology 
for  determining  the  slowness  surface  of  GaAs  was  developed  and  demonstrated.  A  computer  pro¬ 
gram  was  developed  to  predict  the  beam  profile  magnitude  and  phase  as  the  beam  propagates  and 
very  accurate  results  were  obtained.  Future  work  in  this  area  will  focus  on  SAW  propagation 
through  large  electrode  arrays  and  the  continued  development  of  GaAs  SAW  waveguide  structures. 
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Appendix  A 

Solution  Method  for  the  Nonlinear  Differential  Semiconductor  Equations 
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The  mathematical  description  of  the  injection  process  is  given  by  equations  (Al)-(A3).  These 
three  relations  are  Poisson's  equation,  the  current  equation  and  the  continuity  for  electrons,  respec¬ 
tively 


V20  =  -1  Nd+  -  n  +  Nsaw 
€ 

(Al) 

\  =  q/b,nE  +  qDnVn 

( A2) 

1^-  -  I  V  Jn  +  Gn  -  Rn 
cH  q 

(A3) 

where  4>  =  electrostatic  potential 

n=  mobile  electron  density 
=  ionized  donor  impurity  density 
qNSAW  =  acoustic  wave  polarization  charge  density 
Jn  =  electron  current  density 
Gn  =  electron  generation  rate 
Rn  =  electron  recombination  rate 
Dn  =  electron  diffusion  coefficient 
=  electron  mobility 

All  equations  and  terms  involving  holes  have  been  excluded  from  this  formulation  since  the  ACT 
device  is  a  buried  channel  transport  device  and  minority  carriers  can  be  neglected.  A  further 
simplification  results  by  dropping  the  generation-recombination  terms  which  are  not  considered  as 
fundamental  factor  in  the  injection  process.  Hot  electron  effects  will  also  not  be  considered  because 
the  fields  which  act  on  the  carriers  are  much  less  than  the  velocity  saturation  field  for  GaAs.  For 
channels  which  are  much  wider  than  one  acoustic  wavelength  a  uniform  distribution  for  the  poten¬ 
tial  and  carrier  concentration  can  be  assumed  in  the  transverse  direction.  This  allows  for  a  two 
dimensional  analysts  instead  of  a  three  dimensional  analysis,  the  final  simplification  involves  the 
time  variable.  Although  injection  is  a  time  dependent  process,  its  duration  is  much  longer  than  the 
relaxation  time  of  the  mobile  electrons.  Therefore  with  <3n/<3t  set  equal  to  zero,  (Al)-(A3)  can  be 
solved  at  a  series  of  time  instants  with  the  assumption  that  at  each  time  instant  a  steady  state  con¬ 
dition  prevails.  These  simplifications  allow  equations  (Al)-(A3)  to  be  written  in  a  reduced  form  as 
given  in  (A4)-(A6) 


V20  =  -JL  No  —  n  +  Nsaw 
€ 

(A4) 

Jn  =  q/znE  +  qDnVn 

(A5) 

O 

II 

C 

> 

(A6) 

where  V  =  —  x  +  —  y 

3x  ay 

Equations  (A4)-(A6)  along  with  appropriate 
description  of  charge  injection.  The  solution 

boundary  conditions  provide  the  mathematical 
is  obtained  numerically  by  the  finite  difference 

method. 


Additional  manipulations  of  these  equations  are  performed  prior  to  discretization  in  order  to 
insure  numerical  stability.  By  assuming  0  =  0  +  2.  where  0  is  an  initial  potential  and  8  is  a  correc¬ 
tion  term,  (A4)  is  written  as  a  linear  differential  equation  for  8.  This  is  obtained  by  replacing  n  in 
(A4)  by  n.e  kT  an£l  then  expanding  the  exponential  function  in  a  Taylor  series  with  all 

terms  of  second  and  higher  order  in  8  being  neglected.  n*  and  0n  are  the  intrinsic  carrier  concentra¬ 
tion  and  the  electron  quasi  fermi  level,  respectively.  A  procedure  developed  by  Scharfetter  and 
Gummed16'  is  used  to  rewrite  (A5)  before  substituting  it  into  (A6),  This  has  been  shown  to  be 
necessary  for  stability. Lastly,  the  equations  for  8  and  n  are  normalized  into  a  dimensionless 
form.  These  normalized  equations  are  discretized  with  finite  differences  to  obtain  two  systems  of 
equations  which  together  will  approximate  the  potential.  0.  and  the  carrier  concentration,  H.  at  a 
discrete  number  of  grid  points  throughout  the  device.  These  two  systems  are  given  in  matrix  form 
in  (A7)  and  (A8).  The  potential  0  is  calculated  from  8  as  shown  in  (A9). 


AjS  =  b! 

(A7) 

A2ft  =  b2 

(A8) 

0  =  0  +  8 

(A9) 

The  strongly  implicit  numerical  method  of  Herber  Stone^1'  isused  to  solve  (A7)  and  (A8). 
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ABSTRACT 


A  theory  is  presented  for  the  determination  of  time  and  frequency  domain  responses  of  sense 
electrode  structures  in  buried  channel  travelling  wave  charge  transfer  devices  for  real  time  analog 
and  digital  signal  processing.  Two  fundamental  nondestructive  sense  techniques  are  described, 
current  sensing  and  voltage  sensing.  Analytical  and  numerical  results  are  derived  for  the  responses 
of  several  sense  electrode  structures.  The  results  obtained  are  shown  to  be  in  agreement  with 
experimental  data  obtained  from  a  simple  voltage  sensing  structure. 
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1.  Introduction 


The  need  for  wideband,  real  time  analog  and  digital  signal  processing  functions  has  given  rise 
to  a  variety  of  semiconductor  processors  which  may  be  integrated  monolithically  with  passive  and 
active  semiconductor  devices.  A  wide  variety  of  signal  processing  functions  Tiav  be  obtained  via 
the  realization  of  linear  and  nonlinear  operations,  including  gain,  summation,  multiplication,  and 
delay.  The  charge  transfer  device  (CTD)  was  developed  to  allow  the  achievement  of  analog  delay, 
the  most  difficult  of  these  functions. 

The  CTD  employs  a  travelling  electrical  potential  well  for  the  transport  of  charged  carriers 
through  the  device.  The  carrier  transport  mechanism  may  be  realized  in  a  number  of  ways, 
including  clocked  electrical  potential  barriers  and  wells  (employed  in  the  charge  coupled  device 
(CCD),111)  or  the  propagation  of  a  sinusoidal  surface  acoustic  wave  (SAW)  in  a  piezoelectric 
semiconductor  (employed  in  the  buried  channel  travelling  wave  CTD  *  ). 

In  a  surface  channel  CTD.  specifically  the  surface  channel  CCD.13*  carriers  flow  near  a  material 
interface  formed  by  a  semiconductor  and  an  insulating  layer.  Metal  electrodes  placed  above  the 
insulator  form  a  metal-insulator-semiconductor  (MIS)  capacitor.  Carriers  are  transferred  via  the 
coupling  of  the  potential  wells  formed  by  adjacent  capacitors  held  at  time  varying  potentials.  The 
surface  channel  CCD  suffers  from  bandwidth  limitation  due  to  interaction  of  the  signal  charge  with 
traps  at  the  semiconductor  insulator  interface  and  because  diffusion  rather  than  drift  is  the  charge 
transport  mechanism.  Buried  channel  CTD's*4*  overcome  this  limitation  by  transporting  the  carriers 
not  at  an  interface,  but  rather  by  establishing  a  potential  maximum  (for  electron  transport)  within 
the  semiconductor  so  that  carriers  may  be  transported  in  the  bulk  material.  Such  devices  do  not 
suffer  the  limitations  due  to  surface  effects,  and  drift  fields  play  a  role  in  charge  transport  in  buried 
channel  CCD’s,  enabling  higher  speed  operation.  However,  buried  channel  devices  utilizing  non¬ 
destructive  sensing  (MDS)  suffer  a  bandwidth  limitation  due  to  the  finite  spreading  of  flux  lines 
emanating  from  the  signal  charge.  In  the  buried  channel  travelling  wave  CTD.  the  charge  is  carried 
in  svnchronism  with  a  travelling  eiectric  field,  and  diffusion  plays  no  role  in  the  charge  transport 


The  general  problem  of  non-destructive  sensing  consists  of  computing  the  time  and  frequency 
domain  responses  of  a  set  of  sense  electrodes  to  3  known  charge  distribution  in  the  channel.  The 
knowledge  of  these  responses  is  necessary  for  the  successful  design  of  sense  electrode  structures  in 
buried  channel  travelling  wave  CTD's  for  the  realization  of  signal  processing  functions.  This  thesis 
presents  an  analytical  method  for  the  computation  of  these  responses.  Several  NDS  electrode 
structures  are  analyzed.  Analytical  expressions  for  the  responses  in  the  simple  cases  are  developed, 
and  numerical  computation  is  employed  in  the  more  complex  structures.  Results  of  the  developed 
theory  are  compared  with  experimental  data  taken  on  an  .NDS  electrode  structure  in  a  simple 
charge  transfer  device. 
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II  Theory 


A.  Current  Sensing 

The  structure  in  Figure  la  depicts  a  line  charge  of  linear  density  QL  carried  at  a  depth  d  by  a 
travelling  wave  of  wavelength  Aw  and  frequency  fw  beneath  infinitely  long,  perfectly  conducting 
plates  A  and  B  separated  by  a  distance  g,  and  connected  together  by  a  load  of  impedance  Z  The 
magnitude  of  the  travelling  wave  is  assumed  to  be  large  enough  so  that  the  line  charge  is  rigidly 
carried  a  constant  depth  below  the  surface.  The  dielectric  constant  of  the  infinitely  thick  material 
is  assumed  to  be  much  larger  than  that  of  the  medium  above  the  surface.  For  simplicity,  the 
dielectric  constant  of  the  upper  medium  will  be  taken  to  be  zero.  The  assumption  of  zero  dielectric 
constant  is  equivalent  to  requiring  that  no  electric  flux  lines  leave  the  channel,  or  that  no  charge  can 
be  induced  in  the  surface  gap.  The  zero  dielectric  constant  assumption  :s  valid  when  the  dielectric 
constant  of  the  channel  material  is  much  larger  than  that  of  the  medium  above  the  gap.  The 
approximation  is  valid  for  real  travelling  wave  CTD’s  since  the  channel  material  in  those  devices  is 
a  semiconductor  having  a  high  dielectric  constant.  This  assumption  permits  an  analytical 
computation  of  the  induced  charge  density  on  the  electrodes  for  arbitrary  gap  size  g.  The  electrical 
center  of  the  induced  charge  density  perfectly  tracks  the  line  charge,  which  flows  in  the  x  direction 
at  a  constant  velocity  given  by 

v  =  A  f  (1) 

The  bne  charge  velocity  vw  is  assumed  to  be  much  less  than  the  speed  of  light  so  that 
electromagnetic  effects  may  be  ignored.  As  the  line  charge  reaches  the  region  beneath  the  gap.  more 
and  more  flux  lines  terminate  on  plate  B  and  fewer  on  plate  A.  Thereby  an  image  current  is 
established  in  the  surface  electrodes,  which  must  be  carried  by  the  load  Z.  The  impedance  of  the 
load  is  assumed  to  be  small,  so  that  the  channel  potential  is  unperturbed  by  the  voltage  drop  across 
the  load  while  the  current  flows  through  it.  For  simplicity.  Z  is  taken  to  be  zero  so  that  plates  A 
and  3  mav  be  assumed  to  be  at  the  same  potential.  The  time  dependent  current  per  unit  channel 
width  i,(t)  and  its  Fourier  transform  I0(f)  are  the  fundamental  building  blocks  for  the  calculation 
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Fig  1.  AC  models  of  simple  curreni  sensing  structures 

(a)  Structure  with  arbitrary  size  gap  between  sur^tce  plates 

(b)  Structure  with  inhnitesimal  gap. 
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of  non-destructive  sense  responses. 

The  solution  for  finding  i(t)  begins  by  solving  Poisson's  equation 

y2(ft  =  -P±*:Ll 


(2) 


in  the  channel  subject  to  the  boundary  conditions 

(t>  =  0  I  x  I  >  g/2 


—  =  (j  I  x  I  <  g/2 

dy  5 


(3a) 

(3b) 


Since  the  potential  is  known  from  the  solution  of  Eq.  (2),  the  normal  electric  field  and  the  induced 
charge  distribution  on  the  metal  electrodes  may  be  calculated.  The  normal  electrical  displacement 
Dy  at  the  surface  may  be  found  by  taking  the  normal  derivative  of  the  potential  evaluated  at  the 
surface 


Dy  =  -  €  ^ 
dy 


(4) 


ly=0 


Since  the  plates  are  assumed  to  be  perfectly  conducting,  the  induced  surface  charge  density  is 
simply  equal  to  the  displacement 


p,ND(x.p.d)  =  -  € 


d<t> 


dy  iv=o 


(5) 


where  p  is  the  x  coordinate  of  the  line  charge. 


The  total  charge  per  unit  channel  width  on  plate  B  may  be  found  by  integrating  ps  along  the  entire 
plate 

Q(p.d)  =  f  pIND(x  p.d)dx 


(6) 


PlateB 


The  current  per  unit  channel  width  may  be  found  by  taking  the  time  derivative  of  Q(p,d) 


B-9 


Since  the  charge  is  moving  at  constant  depth,  only  the  lateral  position  p  is  a  function  of  time,  and 


i0(t)  may  be  found  by  application  of  the  chain  rule 


i«(t)  = 


dQ(p,d)  dp 


dp  dt 


The  time  derivative  of  position  p  is  the  constant  velocity  of  the  line  charge.  Eqs.  (1)  and  (8)  may 


be  combined 


i0(t)  =  Xwfv 


dQ(p.d) 


dp  lp=x.  r  t 


where  the  origin  of  time  is  taken  to  be  when  the  line  charge  is  directly  beneath  the  center  of  the 


First,  consider  the  special  case  of  the  current  sensing  structure  having  an  infinitesimally  small 


gap.  shown  in  Figure  lb.  The  problem  consists  of  finding  the  fields  induced  by  the  line  charge  QL  a 


distance  d  below  an  infinite  continuous  conducting  plane.  Assume  the  origin  to  be  the  point  where 


plates  A  and  B  meet.  The  potential  may  be  found  by  applying  the  method  of  images.  The  potential 


below  the  metal  is  the  same  for  the  structure  shown  in  Figure  2,  equal  and  opposite  line  charges  of 


magnitude  QL  placed  symmetrically  a  distance  d  about  the  plane  occupied  by  the  conducting  plate. 


which  is  taken  to  be  at  zero  potential.  The  potential  at  any  point  in  the  channel  is  given  by 


-Ql 

0  =  - In 

27 re 


; 

(y-d)2  +  (x— p'_ 


(y+d)  +  (x-p)‘ 


The  induced  charge  density  from  Eq.  (5)  is  then 


p1ND(x.p.d)  = 


-QLd/7T 


(x— p)2  +  d~ 


The  total  charge  on  plate  B  is  found  by  integrating  pIND(x,p,d)  from  x=0  to  x=oo 


d/X  =0.375 


To  simplify  the  solution  of  the  problem,  the  transformation 


w  =  u  +  jv  =  In  Z  — 1  1 1 1 
z  +  g/2 

is  employed  which  maps  the  upper  half  of  the  z-plane  (z  =  x  +  jy)  onto  a  horizontal  strip  ol  width 
7 r  in  the  w-plane.  as  shown  in  Figure  5.  If  the  conducting  planes  A  and  B  are  taken  to  lie  on  the  x 
^  -g/2,  v=0  and  x  g/2,  y=0  axes,  respectively,  in  the  w-plane  there  is  a  continuous  conducting 
plane  at  zero  potential  at  v=0  and  a  continuous  interface  at  v=7T  where  the  boundary  condition  on 
the  potential  is 

f!£=0  08) 

dv 

If  a  line  charge  QL  is  present  in  the  z-plane  at  x=p  and  y=d  the  location  of  the  charge  in  the  w-plane 
is  given  by 


P“f  +Jd 

Wo  =  In  - - - 

P  +  f  +  jd 


(19) 


Now.  the  potential  due  to  the  line  charge  located  at  w0  in  the  w-plane  is  found,  and  transformed 
back  to  the  z-plane  to  obtain  the  potential  in  the  channel.  Due  to  the  electrostatic  shielding 
provided  by  the  conducting  plane  at  v=0.  the  fields  in  the  region  0  <  v  <  7r  induced  by  the  line 
charge  located  at  w0  are  unchanged  by  any  charge  distribution  or  material  structure  which  may  be 
present  below  v=0.  Therefore,  we  are  free  to  introduce  a  line  charge  — Ql  symmetrically  below  the 
conducting  plane.  This  new  line  charge  is  located  at  w0.  Also,  let  an  additional  boundary 
condition  of  zero  normal  potential  gradient  exist  at  v  =  —v.  We  now  have  a  structure,  shown  in 
Figure  6,  which  is  entirely  symmetric  about  the  v  axis,  except  for  the  signs  of  the  line  charges. 
Since  the  line  charges  are  of  equal  magnitude  and  opposite  sign,  the  v=0  plane  remains  an 
equipotential  surface  of  value  zero,  even  in  the  absence  of  the  conducting  plane.  Iherelore,  the 
problem  is  reduced  to  finding  the  fields  due  to  two  line  charges  placed  symmetrically  between 
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Fig.  10  Normalized  frequency  domain  current  sensing  responses  for  structure  with  d/X, 


0.375  and  nonzero  gap  sizes. 
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3  voltage  sense  structure. 


'  ly 


71 


1 


Vv(f)  = 


j2rfCs 


71  + 


1 


2j  sin 


j27TfCs 


77-fL 


fWXv 


Ii(f) 


(32) 


A  more  convenient  model  for  the  voltage  sensing  structure  is  the  Thevenm  equivalent  circuit 
of  the  circuit  of  Figure  13a,  consisting  of  a  voltage  source  VTH(f)  in  series  with  the  structural 
capacitance  (  s  driving  the  external  load  of  impedance  Z.  The  Thevenin  equivalent  is  shown 
Figure  13b.  The  voltage  Vv(f)  across  the  external  load  is  found  by  the  voltage  divider  equation 

71 


in 


Vv(f)  = 


71  + 


V  TH^f) 


j27rfCs 


(33) 


thus  the  Thevenin  voltage  source  has  strength 


71  -F 


VTn(f)  - 


j2trfCs 


71 


Vv(f) 


(34) 


Eq.  (34)  may  be  combined  with  Eq.  (32)  to  yield 


VTn(f)  - 


fwAw  C 


I,(f) 
s 


sin 


7rfL 


f  x 

1  w/x  w 


7rf  L 


t35) 


which  may  be  written  as  a  product  of  terms 


VTH(f)  - 


I,(f)Lv(f)  l 


Csf 


S1  w  A  w 


(36) 


where  Lv(f)  is  a  dimensionless  factor  accounting  for  the  response  to  voltage  sensing.  This  factor 
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may  be  written 


may  be  normalized  to  dimensionless  quantities 


vThA) 

hov(t)  “ 

fA(t)  *  pAT(t) 

(42) 

Qjcs 

✓ — s 

H 

> 

L 

Hov(f)  =  - 

=  — H  j(f)L  (f) 

(43) 

Ql/Csfw 

Therefore,  the  rollofF  in  the  frequency  domain  due  to  the  finite  size  of  the  voltage  sensing  electrode 
is  found  by  evaluating  Lv(f)  at  the  desired  frequency. 

Closed  form  expressions  for  the  voltage  in  the  time  and  frequency  domains  are  possible  in  the 
special  case  of  very  small  gaps  between  ground  planes  and  the  sense  electrode.  Combining  Eqs.  (43) 
and  ( 16)  yields 


where 


where  Pr(t)  is  the  unit  pulse  function  previously  defined. 

As  in  the  case  of  current  sensing,  numerical  results  are  available  for  the  responses  for 
structures  with  arbitrary  gaps.  Since  the  frequency  dependence  due  to  voltage  sensing  is  a  closed 
form  multiplicative  factor  in  the  frequency  response,  numerical  data  for  the  current  sensing 
frequency  response  need  only  be  multiplied  by  the  appropriate  factor  (as  in  Eq.  (43))  to  obtain 
normalized  voltage  sensing  frequency  responses.  The  inverse  Fourier  transform  of  the  frequency 
response  may  be  numerically  computed  to  obtain  the  normalized  voltage  sensing  time  domain 
response. 

C.  Distributed  Charge  Density  in  the  Packet 

An  actual  packet  of  charge  in  the  channel  of  a  buried  channel  travelling  wave  charge  transfer 
device  occupies  a  finite  space  due  to  the  fact  that  the  confining  potential  is  spatially  continuous 
rather  than  abrupt.  A  rolloff  in  NDS  frequency  responses  is  expected  due  to  the  finite  size  of  the 
charge  packet. 

The  charge  packet  density  and  shape  are  found  by  the  solution  of  Laplace  s  equation 

V2^tot  (x.y)  =  o  (47) 

where  0lol(x.y)  is  the  sum  of  the  empty  well  channel  potential  4>c  and  the  potential  due  to  the 
mobile  charge.  4>..  The  mobile  charge  density  p(x,y)  is  thus  given  by 

p(x.y)  =  e  V20e(x.y) 

An  approximate  solution  for  the  mobile  charge  density  has  been  found 


(48) 


where  I0Cf )  and  ps( f)  are  the  Fourier  transforms  of  i0(t)  and  ps(t),  respectively.  Similarly,  voltage 


sensing  responses  are  given  by 


v(t)  =  v0(t)  *  ^£lps( t) 

VL 


V(f)  -  Vt(f)^p,(f) 

VL 


where  v0(t)  and  V0(f)  are  the  time  and  frequency  domain  voltage  sense  responses  for  a  line  charge. 


As  in  the  line  charge  case,  normalized,  dimensionless  responses  for  current  sensing  (h^t)  and  Hj(f)) 


and  voltage  sensing  (hv(t)  and  Hv(f))  in  the  time  and  frequency  domains  may  be  written 


=  77T-  =  ^  »«<t>  * 

VL*w  VL 


H,(f)  =  ^  ~  HJO  P,< 0 


hv(t)  =  h,(t)  *  PAT(t)  fv 


=  fw  h01(t)»  *  PAT(t) 

VL'  **  w 


Hv(f)  =  H,(f)Lv(f)  =  ^lHoi(f)ps(f)Lv(f) 
Vl 


where  hoi(t)  and  Hoi(f)  are  the  normalized  dimensionless  NDS  responses  due  to  a  line  charge. 


The  Fourier  transform  of  Eq.  (55)  cannot  be  expressed  in  closed  form,  even  for  specific  values 


of  4> cl/^o-  The  transform  can  be  obtained  using  an  FFT  routine.  For  the  Special  case  of  0cl/0o  - 


1/2,  the  time  dependent  normalized  charge  density  and  its  Fourier  transform  are  shown  in  Figures 


15  and  16,  respectively.  Since  the  charge  packet  is  localized  to  a  region  only  1/3  wavelength  long 


for  0cl/0o  =  1/2,  the  rolloff  added  by  the  packet  is  not  large  (less  than  1.5  dB  for  frequencies 


below  the  travelling  wave  frequency).  The  rolloff  in  the  frequency  response  will  be  dominated  by 
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channel  travelling  wave  charge  transfer  device,  whose  a.c.  model  is  shown  in  Figure  17.  The 
structure  has  d/Aw  =  0.34.  g/Aw  =  1/4  and  L/Au,  =  1.  The  charge  load  of  the  mobile  packet  was 
assumed  to  be  50%.  The  voltage  on  the  sense  electrode  was  detected  by  a  high-input  impedance 
amplifier  which  drove  a  network  analyzer  configured  to  measure  insertion  loss. 

The  experimental  and  theoretical  frequency  responses  are  shown  in  Figure  18.  Agreement 
between  theory  and  experiment  is  very  good  for  small  and  moderate  values  f/fw.  The  difference 
between  the  theory  and  experimental  result  reaches  an  error  of  about  1.5  dB  at  f/fw  =  1/2.  at  which 
point  the  theory  predicts  the  frequency  response  to  be  approximately  13.5  dB  below  the  d.c. 
response.  The  major  source  of  error  is  the  uncertainty  in  the  exact  value  of  d/Aw,  which  cannot  be 


directly  measured.  Instead,  it  is  derived  from  the  level  at  which  the  device  is  d.c.  biased,  and 


depends  on  the  physical  properties  of  the  semiconductor,  including  the  doping  density  of  the 


channel  material.  An  uncertainty  in  the  value  of  - —  of  ±  20%  would  give  an  error  in  the 

A  ytr 


frequency  response  at  f/fw  =  1/2  of  ±  1.88  dB.  which  would  account  for  the  difference  in  theory 


and  experiment  seen  at  f/fw  =  1/2.  Additional  error  in  the  agreement  between  theory  and 


experiment  would  result  from  error  in  the  estimation  of  normalized  charge  load,  due  to  uncertainty 


in  the  wave  potential  (f>0 . 


The  theory  successfully  predicts  the  salient  features  observed  in  measured  NDS  responses: 


•  Approximately  linear  decibel  rolloff  due  to  the  spreading  of  flux  lines  of  the  charge  in  the 


channel. 


dependence  due  to  the  finite  size  of  the  v 


se  electrode. 


Amplifier 


Network 

Analyzer 


Electrode 


IV.  Conclusion 


A  theory  is  presented  for  the  computation  of  the  time  and  frequency  domain  responses  of 
simple  non-destructive  sensing  electrode  structures  in  buried  channel  travelling  wave  charge 
transfer  devices.  Two  fundamental  sensing  methods  are  analyzed:  current  sensing  and  voltage 
sensing.  Current  sensing  is  the  detection  of  the  current  consisting  of  the  image  of  the  travelling 
charge  packet  flowing  in  the  buried  channel.  Voltage  sensing  is  the  detection  of  voltage  induced  by 
the  image  charge  on  an  isolated  sense  electrode.  Current  sensing  responses  are  computed 
analytically  in  the  zero  gap  case,  and  numerical  solutions  are  presented  for  the  case  of  arbitrary 
gap.  It  is  shown  that  the  voltage  sensing  response  may  be  expressed  as  the  superposition  of  the 
fundamental  current  sensing  responses  in  the  frequency  domain.  A  numerical  analysis  of  the 
dependence  of  the  NDS  frequency  response  on  the  shape  of  the  charge  packet  is  presented,  in  the 
case  cf  a  one-dimensional  charge  distribution  in  the  channel. 

Experimental  results  for  the  voltage  sensing  frequency  response  are  presented  for  a  simple  one 
tap  voltage  sense  structure,  and  good  agreement  between  theory  and  experiment  is  observed. 

The  theory  presented  in  this  thesis  has  been  shown  to  correctly  predict  the  nondestructive 
sensing  response  for  a  simple  buried  channel  travelling  wave  charge  transfer  device.  It  is  expected 
to  be  a  useful  tool  in  the  analysis  and  design  of  signal  processors  employing  the  buried  channel 
travelling  wave  CTD  concept. 
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2.  THE  PROBLEM  AND  ITS  SEPARATION 


The  basic  problem  addressed  by  this  thesis  is  the  determination  of  the  potential  distribution 
for  a  buried  channel  charge  transfer  device  built  with  reverse  biased  Schottky  plates.  This  type  of 
structure  is  in  common  use  in  GaAs  charge  transfer  devices.  In  particular,  various  types  of  struc¬ 
tures  for  a  buried  channel  traveling  wave  charge  transfer  device  were  analyzed.  Of  course,  it  is 
virtually  impossible  to  solve  any  real  physical  problem  exactly:  thus,  approximations  had  to  be 
made.  The  approximations  used,  however,  are  consistent  with  those  in  previous  published  litera¬ 
ture. 

Figure  1  shows  typical  structures  which  occur  in  but  ed  channel  charge  transfer  devices.  To 
solve,  for  these  potentials  without  use  of  an  approximation  would  be  a  difficult  task.  To  simplify 
the  analysis,  two  approximations  are  made  throughout  this  thesis.  One  is  the  depletion  approxima¬ 
tion  and  the  other  is  the  uniform  depletion  depth  approximation.  Both  of  these  approximations  are 
consistently  used  in  previously  published  literature.  It  is  possible  to  incorporate  a  non-uniform 
depletion  depth  into  the  potential  models,  but  the  calculation  of  the  potential  would  become  a 
cumbersome,  iterative  procedure.  The  first  step  of  this  iterative  procedure  is  the  analysis  given  in 
this  thesis.  Future  works  may  contain  models  that  include  non-uniform  depletion  depths. 

Included  in  Fig.  1  are  periodic  structures  and  non-periodic  structures.  Most  of  the  non¬ 
periodic  structures  must  be  solved  with  a  periodic  algorithm  using  a  long  period.  For  the  isolated 
gap  case,  shown  in  Figs,  la  and  lb.  if  one  can  assume  €2=0,  then  a  periodic  algorithm  need  not  be 
used.  In  this  case,  a  conformal  mapping  solution  can  be  used  instead  of  a  Fourier  series  solution. 
The  conformal  mapping  solution  is  also  used  to  select  the  period  required  in  the  periodic  algori;  «m 
to  approximate  a  non-periodic  structure.  One  might  ask  why  one  should  even  bother  with  a  con¬ 
formal  mapping  solution  when  a  Fourier  series  solution  can  be  used  to  solve  virtually  any  arbi¬ 
trary  structure  The  answer  is  simple.  A  conformal  mapping  solution  is  quicker,  more  accurate, 
and  easier  to  use  than  a  Fourier  series  solution.  Thus,  a  conformal  mapping  solution  should  be 
done  whenever  possible. 


C-7 


*  -  * »  v  v-  “ k  *> 


v‘v\-  V 


fi%  _*«W  JV<  -S  kVi, 
1 1*4" 


Typical  structures  encountered  in  buried  channel  charge  transfer  devices 

(a)  Isolated  gap  with  symmetrical  biasing. 

(b)  Isolated  gap  with  unsymmetrical  biasing. 

(c)  Isolated  sense  electrode. 

(d)  Arbitrary  periodic  structure. 


Fig.  1. 


3.  CONFORMAL  MAPPING  SOLUTIONS 


In  general,  solving  boundary  value  problems  that  have  mixed  boundary  conditions  is  not  an 
easy  task.  For  the  case  of  a  single  isolated  gap.  Fig.  4a,  an  accurate  approximation  to  the  potential 
can  be  found  using  conformal  mapping.  The  accuracy  of  the  approximation  depends  on  the  param¬ 
eters  €j,  e2.  a,  d.  Vi  and  V2.  For  e2  =  0,  Fig.  4b,  the  solution  is  exact. 

The  potentials.  V|(x,y)  and  Vn(x,y),  in  the  two  regions  illustrated  in  Fig.  4a  must  satisfy  the 


following  equations. 


V2V!  =  0 


V2V„  =  0 


V,(x.d)  =  0 


V|(x,0)  =  V||(x.O)  =  V,  .  x  <  — a 


V,(x.O)  =  V„(x.O)  =  V2  .  x  <  a 


0V,(x.O)  _  0  V||(x.O) 

6‘  ay  ~ 62  ay 


-a  <  x  <  a  . 


We  see  that  the  potential  is  specified  on  the  plates  and  that  the  electric  flux  must  be  continuous  in 
the  gap,  -a  <  x  <  a.  If  e2  «  €j,  then  Eq.  (3.6)  becomes 

aVi(x.O) 

€i  -  =  0  ,  — a  <  x  <  a  .  (3.7) 

1  ay 

Thus,  al1  the  boundary  conditions  of  Region  I  are  specified.  The  two  region  problem  of  Fig.  4a  has 
been  reduced  to  the  one  region  problem  of  Fig.  4b.  The  necessary  equations  to  solve  for  V(x,y)  are 

V2V  =  0  (3.8) 


V(x,d)  =  0 


'  /  / 


Structure  of  an  isolated  gap. 

(a)  General  case  (two  region  problem). 

(b)  €2  «  €,  such  that  Ey(x.O)  =  0  for  -a  <  x  <  a  (one  region  problem). 
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Conformal  mappings  used  to  solve  for  VA(x.y). 

(a)  z-plane 

(b)  a rplane 

(c)  ci)2-plane 

(d)  3-plane 
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Conformal  mappings  used  to  solve  for  VB(x.y). 
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VjCx.y) 


Ps 

—  nr  —  real((i>3)  +  real(w3) 

-  (d-y) 

el  +  e2 

IT  \ 

where  the  substitution  Es  =  ps/e i  has  been  made. 


This  solution  is  compared  with  the  Fourier  series  solution,  which  will  be  discussed  in  the  next 
section,  and  the  results  are  presented  in  Table  1.  Table  1  shows  that  Eq.  (3.37)  is  valid  for  ratios 
of  a/d  less  than  or  equal  to  1/3. 


The  potential  of  Fig.  4a  for  the  case  Vj  —  V2  and 


a 

d 


1/3  can  now  be  expressed. 


It  is  the 


sum  of  a  linear  potential 


V, 

— -(d-y)  and  Eq.  (3.37)  with  ps  =  — CjVj/d. 
d 


\t  (  'i  2ei  +  e2  V1  /  , 

V|(x.y)  =  — - -(d-y)  - 


ei  +  *2 


e,  +  6, 


y  j  • 

— —{tt — real(wj)  +  real(w3)) 
7 r 


(3.38) 


Thus,  the  potential  due  to  an  isolated  gap  has  been  solved  for  using  the  solution  to  Fig.  7,  i.e., 


Eq.  (3.35).  Equation  (3.35)  turns  out  to  be  quite  useful  for  analyzing  the  effects  of  gaps  in  an 


TABLE  1 

Table  1.  Verification  of  e2  approximation  for  d=l  and  Es=l.  Only  maximum  potentials  are 
specified. 


Fourier 

6i 

a/d 

€2/€i 

Series 

Approx. 

%  Error 

1/2 

1/13 

,432V. 

.423V. 

2.13% 

1/2 

8/13 

.297 

.282 

5.32% 

1/2 

1 

.243 

.228 

6.58% 

1/3 

1/13 

.302 

.296 

2.03% 

1/3 

8/13 

.205 

.197 

4.06% 

1/3 

1 

.166 

.160 

3.75% 

1/4 

1/13 

.231 

.226 

2.21% 

1/4 

8/13 

.155 

.151 

2.65% 

1/4 

1 

.126 

.122 

3.28% 

1/5 

1/13 

.186 

.183 

1.64% 

1/5 

8/13 

.125 

.122 

2.46% 

1/5 

1 

.101 

.099 

2.02% 

1/7 

1/13 

.134 

.132 

1.52% 

1/7 

8/13 

.090 

.088 

2.27% 

1/7 

1 

.073 

.071 

2.82% 

1/10 

1/13 

.095 

.093 

2.15% 

1/10 

8/13 

.063 

.062 

1.61% 

1/10 

1 

.051 

.050 

2.00% 
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otherwise  uniformly  biased  device.  Problems  that  have  more  than  one  gap  can  be  solved  provided 
that  the  voltages  on  the  plates  are  all  the  same.  Equation  (3  35)  can  be  translated  and  superim¬ 
posed  with  the  original  equation  to  simulate  tne  effects  of  multiple  gaps  in  an  otherwise  uniform 
structure.  It  should  be  noted  that  this  doesn’t  yield  an  exact  solution  and  that  superpositions  of 
Eq.  (3.35)  should  not  be  made  if  the  gaps  are  too  close  together.  A  serious  error  could  result  It  is 
recommended  that  gaps  of  width  2a  should  not  be  superimposed  any  closer  than  a  distance  of  3a, 


Equation  (3.35)  is  plotted  in  Figs.  10-14  for  various  values  of  — .  namely.  1/2,  1/3.  1/5.  1/7 

d 


and  1/10.  In  these  plots.  d=l  and  Es  =  1.  Thus,  the  plots  are  normalized.  These  figures  can  be  used 
as  a  quick  reference  for  the  potentials  of  various  structures.  The  potentials  for  structures  like 
those  in  Fig.  la  and  Fig.  lc  can  be  obtained  from  Figs.  10-14  provided  that  the  appropriate  scaling 
factor  is  used.  The  constraints  on  the  types  of  structures  for  which  Figs.  10-14  can  be  used  are 

0  The  surface  voltage  of  the  structure  must  be  uniform  except  for  the  perturbation 

caused  by  plate  separation. 

Either  a/d  <  1/3  or  e2  «  €,. 


2) 


Given  that  these  constraints  are  satisfied.  Figs.  10-14  can  be  used  to  calculate  the  absolute  potential 
deviation  from  the  one-dimensional  potential  given  by  Eq.  (2.1 )  at  a  particular  point. 

The  appropriate  scaling  factor  (SF)  is 


— e,  d 

SF  =  - - —  Es 


e,  +  e-> 


(3.39) 


where  e,  Es  is  the  unperturbed  charge  density  on  the  plate  located  at  y=0.  If  Va  is  the  unperturbed 
voltage  at  y=0  then  using  Eq.  (2.5) 


r  _  v,  -  vy o) 

ts  - - : - 


C3.40) 


Another  wav  to  find  l:s  is  to  integrate  the  space  charge  density  trom  the  channel  minimum  point  to 
the  surface.  Let  vm  denote  the  channel  minimum  noint.  Then 
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\  FOI  P;FR  SERIES  METHOD 


The  previous  section  concentrated  on  the  solution  to  Laplace's  equation  for  isolated  gap  struc¬ 
tures.  It  is  possible  to  use  the  isolated  gap  cas.  to  simulate  multiple  gap  structures:  however,  this 
may  not  always  be  a  good  approximation.  Thus,  a  method  for  solving  Laplace's  equation  with 
multiple  gap  structures  is  needed.  Fourier  analysis  can  be  used  for  such  a  method. 

Figure  15  shows  a  typical  periodic  structure  which  has  three  plates  per  period.  Let  the  sur¬ 
face  potential  at  y=*0  be: 


V(x.O)  =  a0  + 


la, 

n=l 


n7r  .  u 

cos-j— x  +  bn  sin 


ntr 


where  L  is  the  period  of  the  structure.  The  boundary  conditions  that  need  to  be  satisfied  are 


V(x.d)  =  0 

(4.1) 

V(x.-oo)  =  a0 

(4.2) 

V(x.O)  =  ^  if  xeP, 

(4.3) 

av. ,  ev„  cr 

€|— - — (x.O)  =  e2— —  h  xeG 


6y 


ay 


(4.4) 


where  Pj  is  the  domain  of  points  where  the  ’di  surface  plate  lies.  V,  is  the  voltage  on  the  ith  plate, 
and  G  is  the  domain  of  points  that  fall  in  the  gaps.  Let  R]  represent  all  points  (x.y)  between  y=0 
and  y=d.  and  let  the  potential  in  this  region  be  called  V^x.y).  Also,  let  Rn  represent  all  points  such 
that  y<0  and  let  the  potential  in  this  region  be  called  V[(x.y).  An  expression  for  the  potential  in 
regions  R|  and  Rn  can  be  found  by  applying  the  separation  of  variables  to  the  Laplace  equation. 
Fundamental  solutions  from  this  method  are  of  the  form  sinh(ky)sin(kx),  sinh(ky)cos(kx). 

cosh(ky)sin(kx)  and  cosh(ky)cos(kx).  The  values  of  k  are  restricted  to  k  =  n  =  1,2 . by  the 

periodicity  requirement  of  the  potential.  Another  fundamental  solution  of  the  L.aplace  equation  is 
a  linear  potential  of  the  form  Ay  +  B.  A  linear  combination  of  these  fundamental  solutions  can  be 
used  to  obtain  an  expression  for  the  potential  in  regions  Ri  and  Rn-  The  unknown  coefficients  must 
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then  be  evaluated  to  satisfy  the  boundary  conditions.  Solutions  for  the  two  regions  that  satisfy  Fq. 
(4.1)  and  Eq.  (4.2)  are 


(x.y)eR, 


V„(x.y)  =  a0  +  £e 


cos-^-x  +  bn  sin-^-x  .  (x.y)e 


However,  the  boundary  conditions  of  Eq.  (4.3)  and  Eq.  (4.4)  are  not  satisfied.  These  two  boundary 
conditions  at  y=0  are  of  two  different  types.  Mixed  boundary  value  problems  like  this  are  difficult 
to  solve.  In  fact,  there  is  no  known  deterministic  way  to  solve  for  the  Fourier  coefficients.  a0.  aj.... 
and  b,.  b2 .  Normal  Fourier  analysis  can  not  be  used  as  the  surface  potential  is  not  known  every¬ 

where  on  the  surface  y=0.  Thus,  the  best  that  can  be  done  is  to  try  to  find  a  solution  that  is  a  best 
fit  to  the  boundary  conditions. 


McKenna  and  Schryer  have  done  a  study  of  potentials  for  charge  transfer  devices  [3].  In  their 
analysis,  they  used  a  linear  superposition  of  splines  to  represent  the  charge  on  the  plates.  Then, 
using  a  least-square  optimization  routine,  they  picked  the  combination  of  splines  that  minimized 
the  error  in  Eq.  (4.3)  in  a  least-square  manner.  Also,  in  their  analysis.  Eq.  (4.4)  was  automatically 
satisfied  by  the  splines  used.  However,  truncation  of  the  Fourier  series  of  the  splines  can  lead  to 
errors  in  Eq.  (4.4).  It  is  not  expected  that  this  error  would  be  significant,  however.  There  is  little 
doubt  that  McKenna’s  and  Schryer's  method  yields  accurate  results:  however,  the  computation 
necessary  is  quite  lengthy  and  involved.  A  method  is  to  be  presented  that  needs  little  computer 
memory  and  no  special  software,  but  yields  very  accurate  results  as  compared  with  exact  solutions 
known  using  conformal  mapping  techniques. 


Expanding  ps  at  y=0  can.  indeed,  require  many  Fourier  terms.  However,  expanding  the  poten¬ 
tial  at.  y=0  doesn't  require  nearly  as  many  terms  because  the  potential  is  continuous,  but  the  charge 
density  is  not  necessarily  continuous.  Unless  the  electric  fields  must  be  known  exactly  at  the 


definition  of  the  error  to  avoid  erroneous  results. 


If.  on  the  other  hand,  the  number  of  sample  points  equals  the  number  of  coefficients,  then  one 
has  N  constraining  equations  with  N  unknowns  which  can  be  solved  using  Gaussian  elimination. 
The  programming  complexity  of  solving  N  equations  in  N  unknowns  is  far,  far  less  than  optimizing 
N  unknowns  at  N+m  points:  however,  the  computer  memory  is  not  used  in  as  an  efficient  manner 
as  in  an  optimization  routine.  Thus,  there  is  a  trade-off  between  programming  complexity  and 
efficient  hardware  use. 

For  this  thesis.  N  sample  points  were  used  to  solve  for  N  coefficients.  It  was  found  that  com¬ 
puter  memory  space  was  not  the  limiting  factor  in  our  calculations  but  that  numerical  error  was 
the  limiting  factor.  If  the  sample  points  were  picked  too  close  together,  the  matrices  used  would 
become  singular. 

How  one  chooses  the  sample  points  is  crucial.  The  solution  can  explode  if  the  sample  points 
are  not  chosen  carefully.  Instinct  would  lead  us  to  sample  more  densely  in  the  regions  of  most 
rapid  variation.  This  should  be  avoided,  however,  as  the  solution  that  is  assumed  is  band-limited, 
but  the  true  solution  is  not.  At  the  corners  of  the  Schottky  plates  the  potential  varies  very  rapidly. 
This  rapid  variation  cannot  be  simulated  with  a  truncated  series.  Thus,  a  small  zone  should  be  left, 
near  the  corner  of  every  Schottky  plate,  in  which  no  sample  points  are  chosen.  The  best  results 
that  were  achieved  were  with  uniform  sampling  over  [0,L)  except  for  very  small  zones  near  the 
edge  of  every  Schottky  plate,  in  which  no  sampling  was  done. 

If  an  even  structure  is  analyzed,  then  all  of  the  sine  terms  are  zero.  In  addition,  sample  points 
should  be  chosen  over  the  interval  [0.  -t-].  Sampling  over  the  full  interval  may  result  in  degenerate 
equations. 

Using  the  above  procedure,  it  is  possible  to  get  accurate  expressions  for  the  potential  To 
check  the  accuracy  of  the  Fourier  series  method,  the  potential  for  he  structure  of  Fig.  7  was  ca  cu- 
lated  using  conformal  mapping  and  the  Fourier  series  method.  Parameters  used  were  d  =  1  /im.  E  = 
106  V/m  and  various  ratios  of  a/d  where  a  is  the  half-width  ol  gap  and  d  is  the  distance 
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Abstract 


A  modification  to  accommodate  anisotropic  loss  is  made  to  a  method  developed  by  Murray 
and  Ash  for  computing  the  SAW  slowness  surface  of  piezoelectric  substrates.  The  resulting  algo* 
ri'hm  is  used  to  determine  the  slowness  surface  for  the  leaky  or  pseudo-surface  wave  which  pro¬ 
pagates  on  [lOO]-cut  GaAs  and  for  a  cut  five  degrees  off  from  this.  The  data  thus  obtained  are  used 
in  a  diffraction  theory  computer  program  to  predict  the  shape  of  beam  profiles  as  they  propagate  in 
the  substrate  material  and  accurate  predictions  of  beam  shapes  are  obtained. 


Introduction 


In  recent  years  there  has  been  growing  interest  in  the  development  of  monolithic  signal  pro¬ 
cessors  on  GaAs  substrates.  An  example  of  such  a  device  is  the  Acoustic  Charge  Transport  (ACT) 
device1  which  is  very  simply  a  charge-coupled  device  (CCD)  for  which  the  charge  packets  are  car¬ 
ried  by  the  potential  propagating  with  a  surface  acoustic  wave  (SAW)  rather  than  the  potential 
imposed  by  clock  pulses  applied  to  surface  electrodes.  In  order  that  better  ACT  devices  may  be 
developed,  more  knowledge  about  the  SAW  properties  of  GaAs  is  >  eeded.  Heretofore  GaAs  has  not 
been  of  great  interest  as  a  substrate  material  for  SAW  devices  because  it  is  so  weakly  piezoelectric. 
Its  electromechanical  coupling  constant  is  one-sixtieth  of  that  for  LiNb03  and  one-half  of  that  for 
ST-Quartz.  Because  of  this.  SAW  devices  with  reasonable  efficiency,  bandwidth,  and  processing 
gain  are  not  feasible  on  GaAs.  However,  since  GaAs  is  both  a  piezoelectric  crystal  and  a  semicon¬ 
ductor  it  is  an  attractive  substrate  material  for  ACT  devices. 

Were  the  wave  propagating  in  the  [l  10]  direction  on  [l00]-cut  GaAs  a  Rayleigh  wave  as  it  is 
for  YZ  LiNb03  and  ST-Quartz  one  could  use  the  rather  large  existing  body  of  knowledge  describing 
propagation  in  those  materials.  However,  the  wave  of  interest  in  GaAs  is  a  leaky  or  pseudo-surface 
wave  which  is  a  Rayleigh  wave  coupled  to  a  transverse  bulk  shear  as  shown  in  Figure  1.  For  pro¬ 
pagation  precisely  along  the  [110]  direction  the  transverse  shear  component  vanishes  but  as  the  pro¬ 
pagation  vector  deviates  from  this  direction  the  amount  of  energy  radiated  into  the  bulk  shear 
component  increases  dramatically.  An  important  tool  often  used  to  describe  acoustic  propagation 
in  crystalline  materials  is  the  slowness  surface  which  by  convention  is  the  magnitude  of  Ue.  the 
reciprocal  of  the  velocity  vector)  as  a  function  of  its  direction.  Towards  an  economy  of  words  and 
because  the  same  information  is  represented  I  herein  use  the  term  slowness  surface  to  describe  both 
the  velocity  and  the  attenuation  as  functions  of  propagation  direction.  Theoretical  studies  by 
Penunuri2  and  by  Penunuri  and  Lakin.3  utilizing  the  Christoffel  equations,  predicted  the  slowness 
surface  for  zero-degree  GaAs  as  shown  in  Figures  2  and  3.  Obseive  in  Figure  3  how  quickly  the 
attenuation  increases  as  the  propagation  wave  vector  points  away  from  the  [110]  direction.  In  just 
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two  degrees  the  attenuation  due  to  the  bulk  shear  has  increased  two  ord  T  cf  magnitude!  Note, 
however,  that  this  does  not  represent  all  of  the  attenuation.  It  represents  only  the  acoustic  power 
radiated  into  the  bulk  shear  and  does  not  include  the  normal  viscous  losses  of  the  Rayleigh  wave 
One  should  keep  in  mind,  though,  that  the  attenuation  shown  in  Figure  3  is  a  good  indication  of  the 
increasing  amplitude  of  the  bulk  shear  as  the  wave  vector  deviates  from  the  [  1 1 0]  direction.  To  the 
authors  knowledge  there  has  been  no  experimental  evidence  to  verify  or  nullify  the  theoretically 
derived  attenuation  data  and  this  paper  represents  the  first  publication  of  such  data. 

A  key  point,  and  one  which  motivated  this  study  of  zero  and  five-degree  wafers,  is  that  the 
substrate  upon  which  epilayers  are  normally  grown  are  not  precisely  [lOO]-cut  but  are  instead  two 
degrees  off  from  that  as  shown  in  Figure  4.  The  reason  for  using  the  two-degree  wafer  is  unclear. 
Some  sources  say  that  this  is  essential  for  the  growth  of  high  quality  epitaxial  layers4  and  others 
say  that  the  evidence  as  to  whether  or  not  the  surface  normal  being  off  by  two  degrees  makes  any 
difference  is  inconclusive.  Regardless  of  this,  the  custom  persists  and  epitaxial  growth  is  typically 
performed  on  two  degree  wafers  for  GaAs  and  for  other  III- V  compound  semiconductors. 

The  method  used  to  compute  the  slowness  surface  in  this  paper  is  the  technique  described  by 
Murray  and  Ash7  extended  to  deal  with  the  anisotropic  loss  of  the  leaky  wave  on  Ga  \s.  The  wave 
vector  k(0)  is  assumed  to  be  complex  and  the  resulting  Helmholtz  equation  and  subsequent  Green's 
function  can  be  used  to  derive  the  appropriate  angular  spectrum  of  plane  waves  formulation.  Both 
the  extended  technique  for  computing  the  slowness  surface  and  the  diffraction  theory  used  to  pro¬ 
pagate  beam  profiles  result  from  this  assertion  that  k  be  complex. 


n0  =  whole  number  of  wavelengths  between  z t  and  zm  along  the  z-direction 


and  where  we  have  taken  to  be  in  the  principal  value  range,  —v  <  (0/— </>m)  $  v.  We 

have  that  n’(kx)  is  an  integer  incremented  by  ±  1  when  (0/— 0m)  crosses  ±  tt  and  is  zero  when  kx  = 

0. 

As  can  be  seen  from  equation  (10)  an  accurate  procedure  is  needed  to  measure  the  velocity 
along  the  z-direction  so  that  n0  can  be  obtained.  The  method  used  here,  distinct  from  the  one  used 
by  Murray  and  Ash,7  is  to  make  a  longitudinal  scan  as  shown  in  Figure  6  and  to  count  from  the 
phase  response  the  number  of  cycles.  n0,  between  zm  and  zh  Clearly  from  the  longitudinal  scan  one 
can  also  compute  the  wavelength  and  from  the  RF  drive  frequency  obtain  the  velocity  along  the  z- 
direction.  Now  that  n0  is  known  kzr  can  be  computed  according  to  equation  (10). 

To  minimize  computation  time  the  Fourier  transform  specified  in  equation  (1)  is  carried  out 
using  an  FFT  routine  and  as  such  the  number  of  points  in  the  FFT  and  the  spacing  between  sample 
points  in  the  transverse  scan  sets  the  spacing  between  sample  points  in  k-space.  This  in  turn  sets 
the  value  of  kx  associated  with  each  point  in  the  k-space  representation  ol  the  scan  and  with  each 
value  of  kzr(0)  computed.  We  can  now  use  trigonometry  to  compute  the  real  part  of  the  wave  vec¬ 
tor  and  the  angle  9  associated  with  each  sample  in  k-space  as 

k(i)  =  Vkx2(i)  +  kz2r0)  (W 


0(i)  =  tan  1 


(12) 


where 


i  =  k-space  sample  point  index 

WTe  now  have  the  wave  vector  and  a  relation  between  i  and  9  and  so  we  can  obtain  the  velocity 
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Typical  Phase  Response 


Figure  6.  Longitudinal  scans  for  determination  of  v(0  =  0) 


HI.  Experimental  Procedure 


For  these  experiments  split  finger  inter-digital  transducers  (IDT)  were  fabricated  on  zero  and 

o 

five-degree  Cr-doped  Horizontal  Bridgman  wafers  using  approximately  2000  A  of  Al.  The  fingers 
and  spacings  of  the  IDT  were  1  /xm  and  hence  the  nominal  operating  frequency  was  360  MHz.  The 
transducers  were  30  wavelengths  wide  and  100  wavelengths  long. 

The  laser  probe  system  used  to  make  the  longitudinal  and  transverse  scans  has  been  described 
in  a  concurrent  publication6  and  will  not  be  discussed  here  in  detail.  It  is  sufficient  to  say  that  the 
magnitude  and  phase  of  the  surface  displacement  as  a  function  of  position  is  recorded  and  that  typ¬ 
ically  signal-to-noise  ratios  of  50  to  60  dB  are  achieved.  Also  it  should  be  noted  that  for  the  exper¬ 
iments  described  herein  a  2  /xm  diameter  laser  spot  was  used. 

The  first  step  in  the  experiments  was  to  establish  a  good  estimate  of  the  velocity  along  the 
[llO]  direction  and  this  was  done  by  the  following  procedure.  Before  the  longitudinal  scans  shown 
in  Figure  6  were  made,  the  axes  of  the  translation  stage  upon  which  the  device  is  mounted  are 
aligned  under  computer  control  so  that  the  x-axis  is  parallel  to  the  transverse  features  of  the  IDT 
within  0.01  degrees.  A  sequence  of  three  scans  are  made  along  the  z-direction.  The  first  over  a 
range  of  100  /xm  in  0.1  /xm  steps,  the  second  a  much  coarser  scan  over  4300  /xm  in  1  /xm  steps,  fol¬ 
lowed  by  another  fine  increment  scan  over  the  final  100  /xm  of  the  propagation  path  at  0.1  /xm 
increments.  With  the  phase  response  of  the  fine  increment  scans  one  can  establish  a  first  order  esti¬ 
mate  of  the  wavelength  and  by  counting  the  large  number  of  wavelengths  in  the  coarse  middle  scan 
one  can  refine  this  estimate  further.  Knowing  the  RF  drive  frequency  one  can  then  compute  the 
velocity.  v(0=O)  along  the  [  1 10]  direction. 

Once  the  wavelength  has  been  established  transverse  scans  as  shown  in  Figure  5  can  be  made 
and  the  data  required  for  the  slowness  surface  calculation  described  in  Section  II  can  be  obtained. 
For  the  computations  presented  in  this  paper  the  spacing  between  sample  points  was  2  /xm  and  the 
extent  of  the  transverse  scans  went  well  beyond  the  beam  profile  to  ensure  that  all  of  the  propagat¬ 
ing  surface  wave  was  sampled.  If  one  uses  the  beam  provided  by  the  transducer  to  compute  the 
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In  Figure  10  we  present  the  velocity  curve  for  the  five-degree  wafer  along  with  the  data  from 
the  zero-degree  wafer.  Observe  that  along  the  [l  10]  direction  the  velocity  for  the  five-degree  wafer 
is  2852.01  m/sec  compared  with  2861.9  m/sec  measured  for  the  zero-degree  wafer.  The  asym¬ 
metry  in  the  five-degree  velocity  curve  is  most  probably  due  to  a  lack  of  SAW  energy  at  the  higher 
propagation  angles.  Because  of  the  way  the  xylene-rubber  cement  slit  is  created  there  is  not  a  great 
deal  of  control  over  its  width  ana  hence  the  amount  of  energy  in  the  higher  angles.  The  effect  of 
this  decreased  energy  content  is  also  observed  in  the  attenuation  data  shown  in  Figure  11  where  at 
angles  above  five  degrees  the  curve  becomes  noisy.  If  we  assume  that  the  z-directed  velocity  varies 
proportionally  with  the  angle  of  the  cut  we  would  expect  that  the  velocity  for  the  two-degree 
wafers  used  for  epitaxial  growth  would  be  approximately  2858  m/sec. 

In  Figure  11  we  compare  the  attenuation  data  from  the  five-degree  device  with  that  of  the 
zero-degree  device.  Up  to  about  three  degrees  off  from  the  [  1 1 0]  direction  the  two  curves  compare 
quite  favorably  but  beyond  that  point  the  five-degree  curve  undergoes  a  steep  increase  before  it 
becomes  lost  in  the  noise.  This  would  give  some  indicrlion  that  the  leaky  wave  component  is 
stronger  for  the  five-degree  cut  than  for  the  zero-degree  cut  but  that  along  the  [110]  direction  the 
viscous  loss  is  much  greater  than  the  leaky  wave  loss  for  both. 

The  true  test  of  the  measured  slowness  surface  is  whether  or  not  one  can  use  it  to  accurately 
predict  the  propagation  of  a  beam  profile.  Herein  we  have  developed  a  computer  program  to  pro¬ 
pagate  an  acoustic  beam  a  specified  distance  using  the  calculated  slowness  surface.  The  theoretical 
basis  for  this  program  is  the  angular  spectrum  of  plane  waves  technique  described  in  Section  II. 
The  procedure  is  to  take  transverse  laser  probe  scans  at  positions  1  and  2  as  shown  in  Figure  12  and 
to  enter  the  scan  at  position  1  into  the  computer  program  which  implements  the  diffraction  theory. 
The  user  enters  the  distance  in  the  z-direction  between  scan  1  and  scan  2  and  the  program  performs 
an  FFT  on  scan  1  then  propagates  it  the  prescribed  distance  using  the  previously  computed  slowness 
surface  and  finally  performs  an  inverse  FFT  to  get  back  to  the  spatial  domain.  The  computei  pro- 
paga  ed  beam  is  then  compared  with  the  actual  scan  taken  at  position  2. 
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where  <£bo  is  the  zero  bias  barrier  height  and 


E  = 


(3) 


Nd  is  the  donor  concentration  in  the  epitaxial  layer.  VR  is  the  applied  reverse  bias  voltage.  Vbj  is  the 
built-in  potential,  q  is  the  fundamental  electronic  charge,  and  es  is  the  effective  permittivity  of 
GaAs. 

Based  upon  Equation  (l),  an  increase  in  barrier  height  or  a  decrease  in  the  value  of  the 
Richardson  constant  will  lead  to  a  reduction  in  leakage  current.  Assuming  that  the  Schottky  barrier 
height  is  determined  by  the  metal  work  function.  d>m,  and  is  represented  by 

0bo=^m“Xs  (4) 
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5  Liftoff.  A  standard  liftoff  technique  is  used  [14,15], 

6.  Alloy.  The  ohmic  material  is  alloyed  for  five  minutes  in  a  tube  furnace  at  450°C  with  a 


hydrogen  ambient. 

With  the  exception  of  alloying,  the  Schottky  layer  is  formed  following  the  same  procedur* 

B  Modifications 

1.  Polyimide  Rings 

Use  of  polyimide  rings  requires  an  intervening  process  between  the  ohmic  and  Schottky 
metallizations  (see  Figure  1).  In  order  to  determine  the  effect  of  the  pol^  imide  rings  and  maintain  a 
controlled  experiment,  half  of  the  diodes  on  the  wafer  are  fabricated  with  rings  and  half  without. 
Assuming  that  the  polyimide  rings  will  minimize  the  electric  field  peaking  at  the  periphery,  the 
diodes  with  rings  are  expected  to  have  leakage  current  values  approaching  those  predicted  by 
theory. 

For  the  half  of  the  wafer  with  rings,  the  polyimide  processing  is  as  follows.  To  avoid 
pinholing.  two  0.5  fim  thick  layers  of  polyimide  are  spun  on  at  4000  rpm  for  60  seconds.  The 
polyimide.  type  PI  2555.  is  diluted  1:1  with  T9035  thinner.  Bake  times  and  temperatures  are 
optimized  to  produce  ideal  edge  profiles.  In  this  case,  a  45°  slope  is  desired.  The  first  layer  is  baked 
at  150°C  for  20  minutes,  the  second  at  140°C  for  30  minutes.  A  laver  of  AZ4110  positive 
photoresist  is  then  spun  on  and  baked  at  90°C  for  30  minutes.  The  wafer  is  exposed  to  ultraviolet 
light  for  45  seconds,  then  developed  in  a  0.22N  tetramethylammonium  hydroxide  (TMAH) 
solution  that  causes  overcutting  of  the  polyimide  [  16].  As  mentioned  previously,  the  polyimide 
bake  times  and  temperatures  determine  the  amount  of  overcutting. 

The  TMAH  solution  provides  a  one-step  process  for  developing  first  the  photoresist,  then 
etching  the  polvimide.  A  spray  rinse  with  acetone  is  used  to  remove  the  photcresis'.  from  the  r;ngs. 
High-temperature  bakes  (one  hour  at  180°C  then  two  hours  at  225°C)  cure  the  polyimide  and 
render  it  resistant  to  the  effects  of  further  processing. 
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To  eliminate  nonuniformities  between  wafers,  each  of  the  wafers  containing  diodes  with  a 
particular  metallization  is  sawed  into  four  sections.  Therefore,  the  diodes  of  each  metallization 
have  undergone  the  same  processing.  The  sections  contain  samples  of  all  three  diode  areas,  and  each 
section  is  annealed  for  60  minutes  at  temperatures  of  125,  1 75,  225.  or  275°C. 


IV  TESTING 


A.  Test  Procedures 

All  leakage  current  measurements  are  made  with  an  Alessi  Industries  wafer  prober  and  a 
Tektronix  577  curve  tracer.  Since  currents  are  in  the  nanoamp  range,  measurements  must  be  made 


in  the  dark.  Light  shining  on  the  sample  will  cause  optical  generation  of  carriers  which  will 
drastically  increase  the  leakage  current  values. 


In  addition  to  measuring  leakage  currents,  barrier  height  values  are  determined  using 
current-voltage  ( 1— V )  measurements  [22],  Current  density  (Js)  is  plotted  versus  forward  voltage  in 
the  linear  region  of  operation,  and  an  extrapolated  value  of  Js  at  zero  volts  is  obtained.  Barrier 
height  is  then  calculated  with  the  following  equation: 


<f>bo  =  VT 


(6) 


B.  Results  and  Interpretation 
1.  Polyimide  Rings 

As  expected,  the  diodes  with  polyimide  rings  display  much  lower  leakage  current  values  than 
those  without  rings.  Figures  2  and  3  show  the  increase  of  leakage  current  with  area  and  the 
difference  between  leakage  current  values  of  the  chromium-copper/aluminum  diodes  with  and 
without  polyimide  rings.  The  measured  leakage  currents  fall  relatively  close  to  the  theoretical 
values,  and  since  thermionic  emission  theory  does  not  include  edge  effects,  the  data  indicates  that 
the  polyimide  rings  effectively  minimize  excess  leakage  current  at  the  diode  periphery.  The  diodes 
without  rings  have  higher  leakage  currents,  and  the  current  increases  more  rapidly  at  the  outset 
but  then  starts  to  level  off.  If  the  plot  were  extrapolated  out  to  infinite  area,  it  appears  that  the 
two  curves  would  meet  since  edge  effects  would  be  negligible. 

Based  upon  a  comparison  between  diodes  with  polyimide  rings  and  diodes  without,  the 


i 


leakage  current  can  be  separated  into  a  bulk  leakage  plus  generation-recombination  current 
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Figure  2.  Leakage  current  versus  diode  area  for  the  lower  doped  water 
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component  that  scales  with  area  and  an  edge  leakage  component  that  scales  with  circumference.  An 
equation  for  the  reverse  leakage  current  might  look  like 

Ir  =  (lb  +  ls)  +  le  (7) 

where  Ib  is  bulk  leakage.  Ie  is  edge  leakage,  and  Ig  is  generation-recombination  current.  For  the 
diodes  with  polyimide  rings,  the  edge  leakage  term  is  approximately  zero;  therefore,  the  leakage 
current  scales  linearly  with  area.  In  addition,  if  the  applied  reverse  voltage,  VR,  is  small  enough, 
i.e.,  less  than  one  volt,  then  the  generation-recombination  term  is  proportional  to  the  square  root  of 

Vr. 

Experimental  results  are  higher  than  theoretical  predictions  in  all  cases.  This  is  most  likely 
due  to  generation-recombination  current  in  the  depletion  region,  which  is  neglected  in  theoretical 
calculations.  Figure  4  supports  the  existence  of  depletion-layer  generation-recombination  current 
since  the  difference  between  theoretical  and  experimental  leakage  current  values  increases  with 
increasing  reverse  bias. 

The  barrier  heights,  as  determined  by  1-V  measurements,  remain  constant  with  area,  and  the 
values  are  shown  in  Table  1.  The  measured  barrier  heights  provide  an  accurate  estimate  of  the 
magnitude  of  the  leakage  current  for  the  polyimide  ring  diodes.  The  standard  deviation  for  the 
chromium-copper/aluminum  diodes  without  rings  and  the  chromium-copper/aluminum  diodes 
with  rings  on  the  lower  doped  wafer  is  0.01  eV.  For  the  chromium-copper/aiuminum  diodes  with 
rings  on  the  higher  doped  wafer,  the  standard  deviation  is  negligible.  A  deviation  in  barrier  height 
is  responsible  for  scatter  in  the  data  points,  since  a  0.01  eV  change  in  barrier  height  can  alter  the 
leakage  current  by  as  much  as  a  factor  of  1.5. 

2.  Different  Metallizations 

For  both  doping  levels  investigated  (0.4  X  1013  cm'3  and  1.8  X  1015  cm'1),  the  aluminum 
diodes  have  the  largest  leakage  current  values,  and  the  chromium-copper/aluminum  diodes  have 
the  lowest.  This  variance  in  leakage  current  values  is  displayed  in  Figure  5.  Barrier  height 
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measurements  correlate  with  the  observed  leakage  currents.  The  titanium  diodes  exhibit  leakage 
current  values  between  the  other  two  metallizations.  However,  the  titanium  diodes  fabricated  on 
the  lower  doped  wafer  display  leakage  currents  closer  to  those  observed  with  the  aluminum  diodes, 
whereas  the  leakage  currents  of  the  titanium  diodes  on  the  higher  doped  wafer  fall  closer  to  the 
chromium-copper/aluminum  values.  The  discrepancy  may  be  due  to  variations  in  surface 
preparation  and/or  slightly  different  conditions  during  the  titanium  evaporations. 

The  scatter  in  the  data  points  and  variance  from  one  wafer  to  the  next  indicate  that  surface- 
state  effects  are  a  significant  factor  in  determining  leakage  current.  Auger  analysis  has  revealed  a 
correlation  between  high  surface  leakage  and  the  presence  of  a  surface  gallium  oxide  layer  [23]. 
The  mechanisms  proposed  for  the  excess  leakage  are  band  bending  due  to  charges  in  the  surface 
oxide  and  conduction  through  a  thin  layer  of  arsenic  located  between  the  surface  oxide  and  the 
gallium  arsenide  [24], 

3.  Annealing 

Different  trends  are  observed  for  each  metallization  as  annealing  temperature  is  increased. 
Leakage  current  decreases  with  increasing  temperature  for  the  aluminum  diodes  and  decreases  after 
reaching  a  certain  temperature  for  the  chromium-copper/aluminum  and  titanium  diodes. 

The  polyimide  ring  structure  with  the  chroir. ium-copper/aluminum  Schottky  metallization 
yields  less  diodes  suitable  for  controlled  annealing  experiments  since  half  of  the  diodes  have  rings 
and  half  do  not.  Theretoic.  annealing  of  the  chromium-copper/aluminum  diodes  is  only  performed 
at  125°C  and  175°C.  Although  I-V  measurements  indicate  an  increase  in  barrier  height,  leakage 
current  values  are  greater  after  annealing.  Figure  6  shows  this  trend.  Outdiffusion  of  the 
chromium  or  copper  is  a  possible  explanation  for  the  increase  in  leakage  current. 

The  aluminum  diodes,  on  the  contrary,  show  a  marked  decrease  in  leakage  current  alter 
annealing  at  temperatures  in  the  125-275°C  range,  as  can  be  seen  in  Figure  7  1  his  decrease  in 

leakage  current,  along  with  a  corresponding  increase  in  barrier  height,  agrees  with  the  results  of 
previous  researchers  [25,26],  Evidence  exists  that  indicates  that  a  thin  AlxGa]_xAs  region  forms  at 
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Figure  S.  l  eakage  current  versus  applied  reverse  bias  voltage  for  the  titanium  diodes. 


V.  CONCLUSION 


In  this  thesis  methods  of  reducing  gallium  arsenide  Schottky  barrier  diode  reverse  leakage 
current  are  investigated.  Reverse  leakage  current  is  described  as  the  sum  of  a  bulk  leakage  plus 
generation-recombination  current  component  that  is  proportional  to  area  and  an  edge  leakage 
component  that  is  proportional  to  circumference.  The  polyimide  ring  configuration,  with 
chromium-copper/aluminum  as  the  Schottky  metal,  achieves  the  lowest  (0.1-1. 4  nA)  leakage 
current  values  since  the  edge  component  of  leakage  current  is  minimized.  The  aluminum  diodes 
consistently  have  the  largest  amount  of  reverse  leakage,  but  upon  annealing,  a  drastic  reduction  in 
leakage  current  is  observed.  The  titanium  diodes  display  leakage  currents  in  the  nanoamp  range; 
however,  there  is  significant  variation  between  wafers.  Based  upon  this  data,  reducing  edge  effects 
will  result  in  Schottky  barrier  diodes  with  the  lowest  reverse  leakage  current  values  and  allows  a 
correlation  to  be  made  between  those  values  and  the  measured  barrier  height. 


LIST  OF  REFERENCES 


C.  Rhee,  J.  Saltich,  and  R.  Zwernemann,  "Moat-etched  Schottky  Barrier  Diode  D-splaying 
Near  Ideal  I-V  Characteristics,"  Solid  State  Electronics,  Vol.  15,  p.  1181.  1972. 

C.  Rhee  p.  1181. 

J.L.  Saltich  and  L.E.  Clark,  "Use  of  a  Double-diffused  Guard  Ring  to  Obtain  Near  Ideal  I-V 
Characteristics  in  Schottky  Barrier  Diodes,"  Solid  State  Electronics ,  Vol.  13,  p.  857.  1970. 

P.A.  Folkes.  C.C.  Chang,  and  F.  Lane.  "Process-induced  Leakage  Current  in  a  Low  Noise 
GaAs  MESFET."  J.  Electrochem.  Soc.:  Solid  State  Science  and  Technology,  Vol.  132.  No.  6. 
p.  1417.  June.  1985. 

S.M.  Sze.  Physics  of  Semiconductor  Devices,  New  York:  John  Wiley  and  Sons.  Inc.,  p.  258, 
1981. 

N.  Newman.  "Annealing  of  Intimate  Ag.  Al.  and  Au-GaAs  Schottky  Barriers."  Med.  Phys., 
Dec.  1984. 

B.L.  Sharma.  Metal-semiconductor  Schottky  Barrier  Junctions  and  Their  Applications,  New 
York:  Plenum  Press,  p.  4.  1984. 

S.M.  Sze.  p.  850. 

S.M.  Sze.  p.  204. 

J.  Bardeen.  "Surface  States  and  Rectification  at  a  Metal  Semiconductor  Contact,"  Physical 
Review,  Vol.  71.  No.  10.  p.  717.  May  15,  1947. 

J.L.  Saltich.  "Effects  of  Pre-  and  Post-annealing  Treatments  on  Silicon  Schottky  Barrier 
Diodes."  Proc.  of  the  IEEE,  Vol.  58.  No.  3.  p.  492,  March.  1970. 

A.K.  Srivastava  and  B.M.  Arora.  "Effect  of  Annealing  on  the  Richardson  Constant  of  Al- 
GaAs  Schottky  Diodes."  Solid  State  Electronics ,  Vol.  24,  No.  11,  p.  1049,  1981. 


ASO/ ENEGA 

WRIGHT-PATTERSON  A F8  OH  4543  3 


1 


I 


ASD/ AFALC/ AXAE 

WRIGHT-PATTERSON  A FB  OH  4543  3 


1 


ASD/ENAMA  1 

WRIGHT-PATTERSON  A  FB  OH  45430 


ASD/XRS  1 

WRIGHT-PATTERSON  A  FB  OH  45433 


AFIT/LDEE  -  TECHNICAL  LiBRARY  1 

BUILDING  640*  AREA  B 
WRIGHT-PATTERSON  A  FB  OH  45433-6583 


AFWAL/MLPO  1 

ATTN:  DR.  G.  E.  KUHL 

WRIGHT-PATTERSON  A  FB  OH  4543  3-6533 


AFWAL/MLPO  1 

WRIGHT-PATTERSON  A  FB  OH  4543  3-6533 


AFWAL/MLTE  1 

WRIGHT-PATTERSON  A  FB  OH  45433 


AFWAL/FIES/SURVIAC  1 

WRIGHT-PATTERSON  A  FB  OH  4545  3 


DL-4 


IANDING  OFFICER 
L  AVIONICS  CENT 
ARY  -  0/765 
AN  A  POL  I S  IN  462 


IANDING  OFFICER 
L  TRAINING  EQUI 
NICAL  INFORMAT  1 
DING  2068 
N DO  FL  3281  3-71 

ANDE  R 

L  OCEAN  SYSTEMS 
:  TECHNICAL  LI 

DIEGO  CA  92152- 


AVAL  WEAPONS  CE 

:  TECHNICAL  LI 

A  LAKE  CA  9355  5 


SUPERINTENDENT  (CODE  1424) 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93943-5100 


COMMANDING  OFFICER 
NAVAL  RESEARCH  LABORATORY 
CODE  2627 

WASHINGTON  DC  20375 


NAVELEXSYCOM 
PDE-1  10-33 
WASHINGTON  DC  20363 


REDSTONE  SCIENTIFIC  INFORMATION  CENTER 

US  ARMY  MISSILE  COMMAND 

REDSTONE  SCIENTIFIC  INFORMATION  CENTER 

ATTN:  DRSMI-RPRD 

REDSTONE  ARSENAL  AL  35898-5241 

Advisory  Group  on  Electron  Devices 
Attn:  Security  Office 

201  Varick  Street*  11th  Floor 
New  York  NY  10014 


FRANK  J  SEILER  RESEARCH  LAB 
F J  SRL/NHL 

ATTN:  MAJOR  TERRENCE  F.  DEATON 

US  AIR  FORCE  ACADEMY 

COLORADO  SPRINGS  CO  80840*6528 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 
ATTN:  REPORT  LIBRARY 

MAIL  STATION  5030 
LOS  ALAMOS  NM  87545 


AIR  FORCE  ELEMENT  (AFELM) 
THE  RAND  CORP 
1700  MAIN  STREET 
SANTA  MONICA  CA  90406 


A  E 0 C  LIBRARY 
MS-100 

ARNOLD  AES  TN 


AM/TE 
GRAVE 
5635 


HQ  AFCC/DAPL 
BLDG  P-40  NOR 
SCOTT  A F B  IL 


WASH  DC  20591 


1  839  EIG/EIEE 


JTFPMO 

Attn: 


McLean  V A  22102 


T  OF  TRANSPORTATION  LIBRARY 

l 

V, 

Ns  M-493.2#  Room  930 

CE  AVE.  S.W. 

3? 

(KENNETH  W.  IRBY) 

1 

A 

V 

39534-6348 

S-\ 

>' 

Yr 

1 

f>* 

39534-6348 

fcY 

& 

1 

rr 

at  Director 

Research  Drive 

2 

& 

IN 


V  v\  s 


Vv'VW*. 


DIRECTOR 
NS  A/ CSS 
ATTN:  R 1 

FORT  GEORGE  G  ^  E A  D  E  YD  2075  5-6000 


D I  RECTOR 
NS  A/ CSS 
ATTN:  R  31 

FOPT  GEORGE  G  MEADE  MD  20755-6000 


DIRECTOR 
NS  A/ CSS 
ATTN:  R  4 

FORT  GEORGE  G  MEADE  MD  20755-6000 


DIRECTOR 

NS  A/CSS 
ATTN:  R  5 

FORT  GEORGE  G  WEAD  E  MD  20755-6000 

DIRECTOR 
N  S  A  /  C  S  S 
ATTN:  R8 

FORT  GEORGE  G  MEADE  MD  20755-6000 


D I  RECTOR 
NS  A/CSS 
ATTN:  R  9 

FORT  GEORGE  G  MEADE  MD  20755-6000 


director 
NS  A/ CSS 
ATTN:  S031 

FORT  GEORGE  G  MEADE  MD  20755-6000 


DIRECTOR 
NS  A/CS  S 
ATTN:  S21 

FORT  GEORGE  G  MEADE  MD  20755-6000 


DL-11 


LXfLtfA*  rJi  AX  rx  pul  rxrxr**  ax  -x  rj 


DIRECTOR 
NSA/CSS 
ATTN:  V 30  7 

FORT  GEORGE  G  MEADE  MD  20755*6000 


D I  RECTOR 
NSA/CSS 
ATTN :  W07 

FORT  GEORGE  G  MEADE  MD  20755*6000 


DIRECTOR 

NSA/CSS 

ATTN:  W314  (MRS.  JACKSON) 

FORT  GEORGE  G  MEADE  MD  20755*6000 


D I  RECTOR 
NS  A/ CS  S 
ATTN:  R521 

FORT  GEORGE  G  MEADE  MD  20755*6000 


DIRECTOR 
NSA/CSS 
ATTN:  R  52  2 

FORT  GEORGE  G  MEADE  MD  20755*6000 


DIRECTOR 
NSA/CSS 
ATTN:  R523 

FORT  GEORGE  G  MEADE  MD  20755*6000 


DIRECTOR 

NSA/CSS 

ATTN:  R  5  3  (JOHN  C.  DAVIS) 

9800  SAVAGE  ROAD 

FORT  GEORGE  G  MEADE  MD  2G755-6000 


DL-12 


Wit*. 


Director 

USA  Signal  Warfare  Center 
attn:  AMSEL-SW-EA/  D.A.  Harkleroad 
Vint  Hill  Farms  Station 
War  renton/  VA  22186 

AFWAL/AA Al-2 

Attn:  ILt  G.  Mikulanicz) 

Wright-Patterson  AFB  OH  45433 


Signal  Science  Inc. 

(Attn:  John  P.  King) 

7223  Parkway  Dr. 

Suite  105 
Ha  nover /  MD  21 076 

MIT  Lincoln  laboratory 
(Attn:  Irvin  G.  Stiglitz) 

P.0.  Box  73 
Lexington^  MA  02175 


Birch  Associates  Inc. 
(Attn:  J.  N.  Birch) 

1001  Spring  St. 

Suite  20 6 

Silver  Spring*  MD  20910 

Pol  lard  Road  Inc. 

(Attn:  V.  L.  Lynn) 

Suite  708 

2361  Jef f erson-Davi s  Hwy 
Arlington/  V A  22202 

IDA/CRD 

(ATTN:  L.  A.  LIP0RACE) 

Thanet  Rd. 

Princeton/  NJ  08540 


D ARPA/ST0 

(Attn:  Dr.  Ken  Ando) 

1  400  W  i  l  son  £K  vd 
Arlington  VA  22209-2308 


Dr.  George  Simons 
Harry  Diamond  Labs 
2800  Powder  Mi  It  Rd 
(Attn:  DELHD-RT-RB) 

Adel  oh i  *  MD  20783 

US  Army  Night  Vision  and 
Electro-Optical  laboratory 
Attn:  DELNV-AC  Mr.  Terry  Jones 

Advanced  Concepts  Division 
Ft  Belvoir*  V A  22060-5677 

Mr.  Walter  Gelnovatch 

Dioctor*  US  Army  Electronics  Tech. 

and  Devices  Laboratory  (ERADCOM) 

Attn:  DELET-M 

Fort  Monmouth*  NJ  07703 

Dr.  Edwin  T.  Hunter 

Director*  US  Army  Electronics  Tech 

and  Devices  Laboratory  (ERADCOM) 

attn:  DELET-I 

Fort  Monmouth*  NJ  07703 

Mr.  Harper  Whitehouse 
Naval  Ocean  Systems  Center 
Attn:  Code  741 

271  Catalina  Boulevard 
San  Diego*  CA  92152 

Dr.  Timothy  i  ’fee 
Director  of  Research 
Naval  Research  Laboratory 
4555  Overlook  Drive*  SW 
Washington*  DC  20375-5000 

Mr.  Robert  Evans 
Naval  Research  Laboratory 
4555  Overlook  Drive*  SW 
Attn:  Code  5752 

Washington*  DC  20375-5000 

Mr.  Nelson  Head 

Naval  Research  Laboratory 


Mr.  Bob  Graziano 
Technology  Service  Corporation 
Information  Systems  Division 
1945  Gallows  Road/  Suite  404 
Vi  enna/  VA  2  21  80 

Dr.  Soong  Hak  Lee 

Westinghouse  Defense  Electronics 

P .  0 .  Sox  746 

Mail  Stop  1112 

Baltimore/  MD  21203 

Dr.  Harry  Salvo 

Westinghouse  Defense  Electronics 
P.0.  Box  746 
Mail  Stop  11 12 
Baltimore/  MD  21203 

Mr.  Louis  Lome 

Central  Intelligence  Agency 

Attn;  ORD 

Washington/  DC  20505 


0  I  RECTOR 
NS  A/CSS 
ATTN:  L-222 

FORT  GEORGE  G  MEAD  MD  20755-6000 


Genera  l  Dynami c  s 

attn:  Dr  Mathis/  MS-7202-H 

P.0.  Box  85310 

San  Diego/  CA  92138 

D  ARPA/NTO 

attn:  Dr  Al  Ellenthorpe 
1  40  0  Wilson  Blvd 
Arlington/  VA  22209-2308 


DARPA/DSO 

attn:  Mr  Sven  Roosild 
1  40  0  Wilson  Blvd 
Arlington/  VA  22209-2308 


DL-24 


1 


General  Electric  Company 
attn:  Stephen  Wanuga 
Electronics  Lab/  Rm  244 
Electronics  Park 
Syracuse/  NY  1  32  21 

Hughes 

Industrial  Electronics  Group 
attn:  Preston  Campbell 
P.0.  Box  2940 
Torrance  CA  90509-2940 

Loc  kheed 

attn:  Robe  r t  Ward 
3251  Hanover  St. 

Palo  Alto/  CA  94304 


Teledyne  MMIC 
at  t :  Lynwo  od  Cosby 
1501  Wi  l  son  Bl  vd 
Arlington/  V A  22  209-25  50 


1 


L- 


Rensselaer  Polytech  Institute  1 

at  tn  :  Panka j  Oas 
Troy/  NY  12181 


DBA  Systems 

attn:  William  Druyun 

11781  Jackson  Memorial  Highway 

Fairfax/  V A  22033-3360 


Hughes 

Ground  Systems  Group 
attn:  Davi  d  Ef  f  inger 
P.0.  Box  3310 
Fullerton/  CA  92634-3310 


1^7*  V 


■< 


bm 


iSR 

jSeis 


Grumman  Electronics  Systems 
attn:  Howard  Engleson 
Bethpage/  NY  11714-3585 


DL-25 


Syracuse  University 
attn:  Philipp  Kornreich 
103  Redwood  Drive 
Syracuse#  NY  13201 


Hughes/  Microwave  Products  Div 
attn;  Chia  Chuan  Lee 
P.0.  Box  2940 
Torrance/  CA  90509-2940 


GTE/  Advanced  Components  Technology  Center 
attn:  Harr”  Lockwood 
40  Sylvan  Rd 
Waltham.  M  A  02  254 


S  per  r  y 

attn:  Ken  Matys i k 

Sperry  Park/  P.0.  Box  64525 

St.  Paul/  MN  55164-0525 


Northrop/  Defense  Systems  Div 
atttn:  Fred  Me  Kean 
600  Hicks  Rd 

Pollig  Meadows/  IL  60008-1098 


We st inghouse/  Defense  Group 
attn:  D.  Merge  r i an 
Frriendship  Site/  box  1693 
Baltimore/  MD  21203 


Lockheed-  Georgia  Co 
attn:  R.  W.  Milling 
Marietta/  GA  30063 


Raytheon/  research  Div 
attn:  Gary  Monv  res  s 
131  Spr i ng  St 
Lexington/  MA  02173 


DL-2  7 


1 


Van' an  Associates 
attn:  Michael  Wright 
61  1  Hansen  Way 
Mai l  Stop  K-103 
Palo  Alto*  CA  94303-0883 

Texas  Instruments  1 

attn:  C l aude  Tew 

1  3500  North  Central  Expressway 

P. 0.  Box  2  25936 

Dallas*  TX  75265 

Adams-Russel  Electronics  Co  1 

attn:  Thomas  Zeltwanger 

Communications  Program  Div 

1370  Main  St 

Waltham*  M  A  021 54 

RCA  AT  l abs  1 

attn:  Wolf  Altman 
Moorestown  Corporate  Center 
Route  38 

Moorestown*  NJ  08057 

Draper  Lab  1 

attn:  Paul  VanBroekhoven 
555  Technology  Square 
Cambridge*  ,MA  02139 


Fairchild  1 

attn:  Thomas  Brown 

Space  and  Electronics  Comppany 

Germantown*  MD  20874 


Decision-Science  Applications 
attn:  Arthur  BruckHeim 
1901  North  Moore  St 
Suite  1000 
Arlington*  V A  22209 


DL-29 


.  v.-s,’ 


MISSION 


Rome  Air  Development  Cer<  ter 

RA VC  plans  and  executes  research,  development,  tut  and 
A  elected  acquisition  programs  in  support  o f  Command ,  Control 
Communications  and  Intelligence.  (C3 7)  activities.  Technical 
and  engineering  support  Mithin  areas  of  technical  competence 
is  provided  to  ESV  Program  0 Alices  ( POs )  and  other  ESV 
elements.  The  principal  technical  mission  areas  are 
communication electromagnetic  guidance  and  control,  .  u r- 
velllance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  rnicroimve 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


sr  sxsr 


